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ABSTRACT 



We measure the relation between galaxy luminosity and disk circular velocity 
(the Tully-Fisher, or TF, relation), in the g, r, i, and 2;-bands, for a broadly 
selected sample of galaxies from the Sloan Digital Sky Survey, with the goal of 
^ I providing well defined observational constraints for theoretical models of galaxy 

^ I formation. The input sample of 234 galaxies has a roughly flat distribution of 

I absolute magnitudes in the range —18.5 > M,. > —22, and our only morphologi- 

^ I cal selection is an isophotal axis-ratio cut b/a < 0.6 to allow accurate inclination 

O ! corrections. Long-slit spectroscopy from the Calar Alto and MDM observatories 

^ I yields usable Ha rotation curves for 162 galaxies (69%), with a representative 

^ I color and morphology distribution. We define circular velocities Vso by evaluating 

Qh| the rotation curve at the radius containing 80% of the i-band light. Observational 

O I errors, including estimated distance errors due to peculiar velocities, are small 

^ [ compared to the intrinsic scatter of the TF relation. The slope of the forward 

I TF relation steepens from —5.5 ± 0.2mag/log]^ok^s~^ in the (7-band to —6.6 ± 

I 0.2 mag/log^Q kms~^ in the 2;-band. The intrinsic scatter is cr ~ 0.4 mag in all 

^ ' bands, and residuals from either the forward or inverse relations have an approxi- 

' mately Gaussian distribution. We discuss how Malmquist-type biases may affect 
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the observed slope, intercept, and scatter. The scatter is not dominated by rare 
outliers or by any particular class of galaxies, though it drops slightly, to a ~ 0.36 
mag, if we restrict the sample to nearly bulgeless systems. Correlations of TF 
residuals with other galaxy properties are weak: bluer galaxies are significantly 
brighter than average in the ^f-band TF relation but only marginally brighter in 
the z-band; more concentrated (earlier type) galaxies are slightly fainter than av- 
erage; and the TF residual is virtually independent of half-light radius, contrary 
to the trend expected for gravitationally dominant disks. The observed residual 
correlations do not account for most of the intrinsic scatter, implying that this 
scatter is instead driven largely by variations in the ratio of dark to luminous 
matter within the disk galaxy population. 

Subject headings: galaxies: kinematics and dynamics 



Introduction 



The observed correlation between luminosity and disk rotation speed ( iTully &: Fisher 



19771 ) is one of the fundamental empirical clues to the physics of galaxy formation, in par- 
ticular to the relation between dark matter halos and their luminous baryonic components. 
The TuUy-Fisher (hereafter TF) relation has b een widely exploited as a distance indica- 



tor, in studies of the cosmic distance scale (e.g., iTuUy fc Fisherl 119771 : lAaronson et al.lll986 



TuUv fc Pierce 



Willick 



1990 



2000 



Freedman et al.l 120011) and the large scale pecu l iar velocity field (e.g . , 



Mathewson. Ford, fc Buchhornlll992l : IWilhck et al.lll997l : ICourteau et al.ll2000r ) . 



The ambitious surveys constructed for such studies have usually focused on a relatively nar- 
row range of galaxy types, typically undisturbed late-type spirals, with the goal of obtaining 
a tight relation that can yield precise distances. However, the small scatter and the mea- 
sured parame ters of the TF relation are adopted as key constraints on galaxy formation 
theor i es (e.g., 



Kauffmann. White, fc Guiderdoni 



1998: Somerville fc Primack 



1999 



1993 



Cole et al.lll994J: IMo. Mao, fc White 



Navarro fc Steinmetzl " 20001 ) . even though these may not 
yet have the detail to predict the precise types of model galaxies. The goal of this paper 
is to measure the TF relation for a broadly selected sample of galaxies, with a focus on 
quantifying rather than minimizing the intrinsic scatter and on measuring the correlation of 
TF residuals with other galaxy properties. 



Our sample of galaxies is drawn from the Sloan Digital Sky Survey (SDSS, lYork et al. 



20001 ). an imaging and spectroscopic survey of the North Galactic Cap and selected regions of 
the South Galactic Cap. The SDSS galaxy spectra are obtained through a 3" diameter fiber, 
so they do not yield reliable estimates of rotation velocities for galaxies that are spatially 
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well resolved. We have therefore obtained long-slit Ha rotation curve data for a sample of 
234 SDSS galaxies using the Calar Alto 3.5-m telescope and the MDM 2.4-m telescope. The 
Calar Alto observations of 189 galaxies constituted a substantial portion of the Calar Alto 
Key Project contribution to the SDSS. Selection from the SDSS redshift survey allows us 
to define our sample based on absolute magnitude rather than apparent magnitude, while 
working in a redshift range 5000 km s~^ ^ v < 15000 km s~^ where peculiar velocities add 
little uncertainty to individual galaxy distances. Equally important, the SDSS provides high 
quality, 5-band imaging for all of the program galaxies, and the connection to the much larger 
SDSS database allows statistical application of our results to, e.g., infer the distributions of 
galaxy potential well depths or angular momenta in the local universe. 

As discussed in §2, our selection criteria produce a roughly flat distribution in r-band 
absolute magnitude over the range -18.5 > Mr > -220 We select galaxies with isopho- 
tal axis ratios b/a < 0.6 so that we can make accurate inclination corrections to rotation 
speeds, but we impose no other morphological selection criteria. We obtain usable Ha ro- 
tation curves for 170 of our 234 selected targets (73%). The axis ratio cut slightly reduces 
the representation of early-type galaxies, and the early-type galaxies that pass this cut are 
slightly less likely to yield usable Ha rotation curves. Nonetheless, relative to most previous 
large TF samples, our sample is more representative of the range of disk galaxy types. In 
addition, the redshift range and accurate photometry make our typical observational errors 
smaller than the intrinsic TF scatter, which is essential for an accurate estimate of this scat- 
ter and a full understanding of residual correlations. In brief, this is a TF sample designed 
for studies of galaxy formation, not for measurements of the Hubble constant or peculiar 
velocities. We presented and discu ssed the scaling relations for a disk-dominated subset of 
this sample in (jPizagno et al.ll2005l . hereafter P05). The selection of inclined galaxies ensures 
a small uncertainty on the inclination corrected velocity width. 

There are a variety of ways to define a galaxy's rotation speed from optical or 21cm 
data, and in general they yield sir ailar but not ide ntical results for the correlation between 
luminosity and rotation speed (see IVerheijenI (l200ll ) for a careful investigation of this issue). 
As discussed in §4, our measure of disk rotation speed is based on the value of an arc- 
tangent fit to the rotation curve (ICourteaul 119971 . h ereafter C97) e v aluate d at a position 
containing 80% of the total i-band flux. The work of iTuUy &: Fisherl (119771 ) used 21cm line 
widths rather than optical rotation curves, so strictly speaking our analysis is not the "TuUy- 
Fisher relation" , but we will follow common practice in using this term for the more general 
correlation between luminosity and gas rotation speed. 



^Throughout the paper, we adopt h = / 100 kms ^ Mpc ^ = 0.7 and quote absolute magnitudes for 
this value. One should add 51og(/i/0.7) for other values of h. All logarithms are base 10. 
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The observational study most similar to ou r own, in terms of a broad sample selection, 
is that of lKannappan. Fabricant. &: Franxl (120021 ). hereafter K02, who measured TP relations 
and residual correlations inU, B, and i?-bands for a sar nple of 68 galaxies w ith morphological 
types Sa-Sd from the Nearby Field Galaxy Sample (jJansen et al.l l2000l ). Our sample is 
larger by a factor of ~ 2.5, and the luminosity distribution is rather different: its absolute 
magnitude histogram is roughly flat in the range —18.5 > Mr > —22, with somewhat higher 
representation on the bright end, while K02's histogram peaks at Mb = —18 and declines 
at brighter magnitudes. The selection of our sample at f > 5000 kms~^ and the use of 
SDSS surface photometry to determine disk inclinations makes our observational errors much 
smaller than those of K02, giving us a better handle on the TP intrinsic scatter. However, 
unlike K02 we do not have literature HI data for many galaxies in our sample, nor do we 
extend our sample to faint luminosities. In comparisons to the voluminous observational TF 
literature, we will mostly focus on K02 and on the work of C97, who analyzed large samples 
of disk galaxy optical rotation curves using methods similar to those adopted here, and 
Verheijen (2001, hereafter VOl), who carried out a detailed investigation of the TF relation 
in the Ursa Major cluster using resolved HI rotation curves and optical and infrared surface 
photometry. 

Numerous papers have used TF results to derive conclusions about the physics of galaxy 



of disk galaxies (e.g. 


. Fabei 


1982 


Gunn 1982: Persic & Sahicci 1988: ICole & Kaisei 


1989: 


Whit 


e & Frenk 1991; 


Kauffmann. White. & Guiderdonil 1993; Cole et al. 
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3isenstein & Loeb 


1996; 


Persic. Salucci. 


& Stel 


1996; ] 


Dalcanton. SDereel. & SummersI 199?!; 


Mo. Mao. & Whitel 


1998; 


Courteau & Ri 


2 1999I: 


Firmani & Avila-Reese 2000l: Navarro & Steinmetz 


2000l:l 
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2007; 


Gnedin et al. 2 


006). Our underlying goal is to provide better observational inputs for 



these kinds of theoretical and modeling efforts. We have therefore paid particular attention 
to characterizing our sample selection and observed quantities as straightforwardly as we 
can, and to estimating the intrinsic scatter in the different SDSS bands. We have already 
presented the scaling relations and inferred baryonic ni ass fractions of a disk-dominated sub- 
set of our full sample in P05, and iGnedin et al.l (120061 ) have used these results to constrain 
the distributions of disk-to-halo mass ratios and spin parameters of spiral galaxies. 



2. SDSS Observations and Sample Selection 



The SPSS use s a mosaic CCD camera (iGunn et al.lll998l ) on a ded icated 2.5-m telescope 
(iGunn et al.l 120061 ) to image the sky in five photometric band passes (jPukugita et al.lll996l ) 



- 5 - 



denoted m, g, r, i, The imagi ng data are redu ced by a series of automated pipelines that 



2001 



2001 



Smith et al 



2002 



perform astrometric cahbrat ion (jPier et aLll2003l ). photome tric data reduction (ILupton et al. 



Ivezic. Lupton. fc Schlegelll2004j : iTucker et al 



Stoughton et al.ll2002), and photometric calibration JPukugita et al. 1996 : Hogg et al. 



20061 ) . We choose our 



TF sample from the main galaxy spectroscop ic sample, wh i ch is selected from the SDSS 
imaging data using the algorithm described by lStrauss et al.l (120021 ). 



As discussed in the introduction, our goal is to investigate a galaxy sample that is as 
representative as possible of the full galaxy population, while keeping our sample complete- 
ness high and our observational errors significantly below the expected intrinsic scatter of the 
TF relation. We have chosen to focus on the absolute magnitude range —18.5 > Mj. > —22. 
Brighter galaxies tend to be early-type systems for which it is difficult to obtain Ha rotation 
curves, and the behavior of the TF relation at fainter magnitudes, while an interesting prob- 
lem in itself, brings in additional complications because of the morphological irregularity and 
low rotation speeds of the galaxies. 

We correct our observed velocities for inclination using observed disk axis ratios, but 
select galaxies in such a way that any possible intrinsic disk ellipticities contribute a negligible 
amount of uncertainty to the correction. The intrinsic ellipticities of disks, thought to b e 
~ 5 — 10% rms depending on galaxy type and band (jZaritsky &: Rixl Il995l : iRydenI |2005| ) . 
are an irreducible source of uncertainty in TF studies without 2D velocity fields, since they 
prevent one from perfectly measuring a galaxy's inclination and thus inferring its deprojected 
rotation speed. We have selected galaxies to have a measured axis ratio b/a < 0.6, so that 
a 5% intrinsic ellipticity would change the inclination corrected rotation velocity V^bs/ sini 
oc K5bs(l — &^/a^)~^/^ by < 5%. Random ellipticities of this order would then add ~ 0.15 
magnitudes of scatter to the forward TF relation for a typical slope L ~ V^^, which is 
smaller than typical estimates of the intrinsic scatter by a factor ~ 2 — 3. Note, however, 
that disk ellipticities also contribute scatter via non-circular motions; see the discussion in 
§5.3 below. Specifically, we base our sample selection on the SDSS r-band isophotal axis 
ratio (isoB_r/isoA_r), which is measured at an isophote of 25 magarcsec"^ where the disk 
typically dominates over the bulge. We ultimately base our inclination corrections on the 
results of our 2-D bulge-disk decompositions described in §4.3 below. The intrinsic ellipticity 
of a disk is not the same as the intrinsic axis ratio, where the intrinsic axis ratio is the observed 
axis ratio when a galaxy is viewed edge-on (see §4.5 equation 3). 



Fukugita et all (jl996[ ) actually define a slightly different system, denoted u' , g' 



i', z', but SDSS 

magnitudes are now referred to the native filter system of the 2.5-m survey telescope, for which the bandpass 
notation is unprimed. 
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We assume tha t typical galaxy peculiar velocities relative to the local large scale flow, 
~ 200— 300 km s~^ fIStrauss &: Willickl Il995l ). introduce only a small uncertainty in galaxy 
luminosities. We therefore select our sample to have galaxy redshifts cz > 5000 km s~^, so 
that the corresponding absolute magnitude uncertainties are AM < 51og(l + 300/5000) = 
0.13 mag. At the same time, we want galaxies to be at least several arc-seconds across, so 
that we can get good morphological measurements from the images and ~2" seeing does not 
seriously degrade rotation curve measurements. Since galaxy intrinsic sizes correlate with 
luminosity, we therefore impose a maximum redshift that depends on absolute magnitude: 
9000 kms"^ for -18.5 > Mr > -20, ll,000kms-^ for -20 > Mr > -21, and 15,000kms-^ 
for Mr < -2li These cuts ensure that all galaxies have a half-light diameter 2r^Q > A"; 
the median value in the sample is 2r5o = 19.5", and the 10% and 90% values are 11.3" 
and 34.8" (based on the i-band bulge-disk decomposition described in §4.3 below). The 
faintest apparent magnitude allowed by these cuts (for Mr = —18 galaxies at 9000 kms~^) 
is r = 17.5, brighter than the SDSS main galaxy spectroscopic limit r=17.7 mag. The median 
apparent magnitude of the sample is r = 15.06 mag. We compute galaxy luminosities using 
SDSS Petrosian fluxes and colors using SDSS model colors, both i^"-corrected to redshift 
z = using Blanton et al.'s (l2003al ) kcorrect_v3. lb. We compute distances u sing the SDSS 



heliocentric redshifts corrected to the rest frame of the Local Group barycenter (IWillick et al. 



19971 ). assuming a cosmological model with ^2^=0. 3, fiA=0.7, and h = 0.7. We incorporate 



a distance uncertainty corresponding to 300 km s ^ when calculating disk scale length and 
l uminosity uncertainties, to account for the typical amplitude of small scale peculiar velocities 
jstrauss fc WiUicklll995h . 

This combination of absolute magnitude and redshift limits gives a distribution of can- 
didates that is roughly flat in absolute magnitude over the range —18.5 > M,. > —22, since 
we search for rarer, brighter galaxies over a larger volume. We do not impose any explicit 
cut at Mr = —22 mag, but we have relatively few brighter galaxies in our sample because 
of the declining luminosity function. Ideally, we would search for all galaxies in a relatively 
narrow shell, e.g. 5000 km s~^ < cz < 7000 km s~^, so that they would be as well resolved 
as possible, and we would use weighted random sampling to obtain a flat distribution in 
absolute magnitude. However, we began our observations in June 2001 when the sky area 
covered by the SDSS was still relatively small, and we needed to go to larger distances for 
brighter galaxies to obtain a sufficient number of targets. We have kept our selection criteria 
fixed throughout the course of the observing program, rather than lower the outer redshift 
limits as the SDSS sky coverage increased. All of the galaxies in our sample are included 



•^These cuts are not exact because of small changes in the SDSS photometry and our selection criteria 
over the course of the survey, but they are a close approximation to our redshift limits. 
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in the SDSS Second Data Release (DR2; lAbazajian et al.ll2004j ). and when possible we have 



taken images and photometric parameters for the final analysis in this paper from the public 
database. The selection of candidates was based on the state of reduction of the imaging 
data available at the time of our observations, so there might be small differences between 
the parameters used for selection and the parameters used for analysis. The only SDSS 
spectroscopic parameter that we use is the redshift, to determine galaxy distances. 

Figure [T] shows the r-band absolute magnitude distribution of our sample. The open 
histogram shows all of the galaxies that we spectroscopically targeted, and the filled his- 
togram shows those galaxies that yielded Ha rotation curves good enough for use in our TF 
study. As discussed in §4.1 below, we characterize galaxies with unusable rotation curves 
as fiag-4, those with rotation curves still rising at the outermost data point as fiag-3, those 
with rotation curves just reaching the turnover region as fiag-2, and those with extended 
fiat portions of the rotation curve as fiag-1. Most fiag-4 galaxies simply lack a sufficient 
amount of extended Ha emission, but in four cases they have an extended emission pattern 
that has too low signal-to-noise Ha emission to allow sensible definition of a rotation veloc- 
ity. Overall, we obtained usable Ha rotation curves (fiag-1, 2, 3) for 170 of our 234 target 
galaxies, or 73%. After the bulge-to-disk decomposition, discussed in §4.3, we compared our 
disk position angle to the SDSS isophotal position angle and removed six galaxies for which 
the difference between the SDSS and our position angles would cause a >10% difference 
in rotation velocity. We also removed two galaxies that have disk axis ratios greater than 
0.7, as derived from our bulge-disk decomposition. The final sample size is 162 galaxies. 
The filled histogram is approximately fiat in absolute magnitude, but there are fewer low 
luminosity systems, in part because we concentrated our early observations on the range 
Mr < —19.5. In our TF analysis, we will quote results both for the sample "as is" and 
for the sample weighted to approximate a truly fiat absolute magnitude distribution over 
the range —18.5 > > —22, since the latter provides a well defined target for theoretical 
models to predict. In practice, this reweighting makes little difference to the results (see 
§5.2). 

Ideally, we would like to obtain rotation velocities for a random subset of all galaxies 
in the absolute magnitude range —18.5 > Mj. > —22. Our only explicit morphological 
pre-selection is the axis ratio requirement b/a < 0.6, which is necessary to allow accurate 
inclination corrections to the measured velocities. Since disks are flatter than bulges, this 
axis ratio cut tends to suppress the representation of early type galaxies in the sample. The 
requirement of obtaining a usable Ha rotation curve imposes another, less well-controlled 
selection on galaxy type — the galaxies in our flnal sample are those with a sufficient amount 
of extended, ionized gas. In principle, we could obtain rotation curves for other galaxies using 
stellar absorption lines, but these require considerably longer exposures, and we did not have 
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enough observing time to get a large sample and go deep enough for absorption line rotation 
curves. Completing the sample with absorption line rotation curves would be a valuable 
follow up to the present study, requiring a comparable amount of observing time. 

Figure [2] shows the distributions of our sample galaxies in the color- magnitude plane, 
Mr VS. g-r. Contours show the l uminosity-weigh t ed col or magnitude distribution for the 
full DR2 sample of galaxies from iBlanton et al.l (j2003bl ). with contours containing 25%, 
50%, and 75% of the DR2 total luminosity density. The four data point types indicate the 
quality of the rotation curve obtained for the galaxies. Figure [2] shows the distribution of 
our galaxies relative to a sample without velocity measurements. Galaxies without extended 
Ha emission are concentrated at low luminosity, where the surface brightness is low, and 
at high-luminosity and red colors, where galaxies are more likely to be gas poor. Galaxies 
with rising rotation curves are also more concentrated at low luminosity. However, there are 
usable data in all regions of the color-ma gnitude plane. In accord with the usual trends for 
luminosity-environment correlations (e.g. iHogg et al.ll2004l ). the high luminosity galaxies are 
predominantly in over-dense regions, while the lower luminosity galaxies are found in less 
dense regions. 

Figure [3] compares the prop erties of our sample to a "contr ol" sample with the same 
Mr distribution from SDSS DR4 (lAdelman-McCarthy et al.ll2006l ). We have selected a large 
set of galaxies from DR4j^ that matches our full sample's Mr distribution (open histogram 
of Figure [T]) but has 22 times more galaxies in each 0.5 magnitude bin, thus minimizing 
statistical fluctuations in the control sample. The cumulative g-r distribution of the full 
DR4 sample is shown as a solid line in Figure [3^. The dotted curve shows the distribution 
after restricting this sample to galaxies with i-band isophotal axis ratio h/a < 0.6; the impact 
of the axis ratio cut is tiny. The short-dashed and long-dashed lines show our flag-1,2,3,4 
and flag-1,2,3 distributions. Comparing the short-dashed and long-dashed lines, we find that 
selecting galaxies according tO Ha detectability does not alter the color distribution. There 
is a small but statistically significant difference between our sample of (flag-1,2,3,4) and the 
DR4 sample with h/a < 0.6. We have been unable to identify the cause of this difference, 
since in principle the two samples should differ only in size, but it is small in any case. 

Figure Wp shows the cumulative distributions of the i-band concentration index (cj = 
rgo/rso), where rgo and contain 90% and 5 0% of the Petrosian flux. A typical "early- 
type" division is q > 2.6 (IStrateva et al.l 120011 ): see Figure [TOl below. Comparing the solid 
and dotted lines in Figure Wp shows that the early-type fraction is reduced when the axis 
ratio cut is applied, but the effect is small. Comparing the short-dashed to the long-dashed 



''Using the Skyserver website for DR4: 



http : //skyserver2 . f nal . gov/ dr4/en/ 
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lines shows that selecting galaxies according to Ha detectability reduces the early-type frac- 
tion, but again the effect is small. We conclude that our axis ratio cut and requirement of 
Ha detectability do not strongly alter the color or concentration distribution of the galaxy 
population when compared to a random sample of galaxies with the absolute magnitude 
distribution shown in Figure [H In this range, the "true" early- types (i.e., ellipticals 
and bulge-dominated SOs) that would be strongly suppressed by the axis ratio cut are rare, 
and the high completeness of the Ha observations leaves little room for further selection. 
Classifying the flag-1,2,3,4 galaxies according to Hubble type, shows that the flag-4 galaxies 
are predominantly Sa with a few, if any, true SO galaxies. 



3. Spectroscopic Observations and Data Reduction 

Spectroscopic observations were carried out at the Calar Alto Observatory using the 
TWIN spectrograph mounted on the 3.5-m telescope, and at the MDM Observatory using 
the CCDS spectrograph mounted on the 2.4-m Hiltner telescope. Observations at Calar Alto 
were carried out between June 2001 and October 2002; 30 nights were allocated in total, but 
many were clouded out. Observations at MDM were carried out between February 2003 and 
April 2004, on a total of 18 nights. A total of 237 spectra were taken, with 52 spectra from 
MDM and 185 from Calar Alto. Spectra for 3 galaxies were repeated at Calar Alto and MDM 
to allow comparison between the two telescopes (see the end of this section). As discussed 
above and in §4.1, 170 of the 234 un-repeated spectra yielded usable Ha rotation curves, 
and eight galaxies were eliminated from the sample based on position angle misalignments 
or axis ratios. 

The initial Calar Alto observations were carried out with total exposure times of 1800 
seconds, chosen to provide a high signal-to-noise ratio on the Ha line. The observing strategy 
was updated in later observing runs to 1200 seconds for bright galaxies, 1800 seconds for 
medium galaxies, and two 1200 second exposures for faint galaxies. Most MDM observations 
took three 1200 second exposures. The last MDM run, on which we obtained spectra for 18 
galaxies, took an initial exposure for 1200 seconds. Galaxies that showed no Ha in the first 
exposure were not re-observed. If Ha was detected in the initial exposure, then three more 
1200 second exposures were taken, making the total exposure time 4800 seconds. Table [1] 
summarizes the main characteristics of both spectrographs and the spectrograph setup. The 
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data were reduced using standard IRAI^ routines and the XVistal^ package. 

The raw MDM and Calar Alto data were reduced in a similar manner. The data were 
corrected for bias and dark currents and flat fielded following procedures outlined in the Kitt 
Peak Low-to-Moderate Resolution Optical Spectroscopy Manual0 Wavelength calibration 
and linearization were performed using neon-arc lamps. Cosmic rays were removed from the 
MDM data through median filtering of multiple images and the IRAF routine .CRREJECT., 
or by hand for the Calar Alto data. The telluric lines were also used to test the alignment 
and wavelength calibration of the spectra. 

Accurate measurements of velocity centroids at each spatial position along the slit re- 
quire accurate removal of the galaxy continuum. A noise image was made that included the 
galaxy continuum to accurately account for the Poisson noise due to continuum plus Ha 
photons. The galaxy continuum was removed by averaging ~10A of data on each side of 
the emission lines, and subtracting this from the area containing the Ha emission. Telluric 
emission lines were subtracted from the data by averaging rows in the spatial direction, 
where there was no galaxy data, and subtracting that from the area of th e CCD containing 



the ga laxy data. Tests using the wavelengths of telluric lines tabulated by lOsterbrock et al. 



(119961 ) show that the dispersion axis was aligned to be perpendicular to the columns of the 
CCD to an accuracy of 0.1 A. Four flat-fielded, wavelength calibrated, and linearized spectra 
are shown in Figure HJ illustrating varying levels of Ha detectability and spatial extent. 

Figure O compares the rotation curves for the galaxy SDSS J024459. 89-1-010318. 5, ob- 
served at both Calar Alto and MDM. The overall shapes of the rotation curves are similar, 
with slightly lower velocities in the MDM data. The two arc-tangent function fits (see §4.1 
below) are similar. The arc-tangent fit velocities at 18.4" along the arc-tangent functions, 
the radius containing 80% of the i-band flux (see §4.2), differ by 4.7 kms~^, about 1.3a 
(MDM= 186.0 ± 3.3 kms'\ Calar Alto= 190.7 ± 3.5 kms"^). The other two galaxies with 
MDM and Calar Alto spectra show even better agreement in measured rotation speeds than 
the example shown in Figure [51 We used the Calar Alto rotation curve for this galaxy, 
because the systematically lower MDM rotation curve suggests a modest slit mis-alignment. 
Regardless of these small differences there is still overall good agreement between the Calar 
Alto and MDM data. 



^IRAF is written and supported by the IRAF programming group at the National Optical Astronomy 
Observatories (NO AO) in Tucson, Arizona. NOAO is operated by the Association of Universities for Research 
in Astronomy (AURA), Inc. under cooperative agreement with the National Science Foundation. 



^XVista can be found at: http://ganymede.nmsu.edu/holtz/xvista 



^The manual can be found at: 



http://www.noao.edu/kpno/manuals/12mspect/spectroscopy.html 
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4. Measuring Rotation Velocities and Inclination-Corrected Luminosities 

4.1. Rotation curve fitting 

The procedure to extract rotation curves from tlie spectra is similar to that used by 
C97. The 2-D spectra for a galaxy are extracted into 1-D linear spectra, with each 1-D 
linear spectrum along the spatial direction of the slit being 1 pixel wide (0.41" for MDM and 
0.55" for Calar Alto). To avoid assuming an implicit shape of the emission lines, we measure 
the intensity weighted velocity centroid of the Ha line in each 1-D spectrum. The uncertainty 
in the Ha centroid is measured using the signal-to-noise ratio (SNR) at each pixel with Ha 
flux. The Ha line centroid has typical uncertainties of 2 — 12 kms~^, depending on the total 
SNR of the emission line. 

Following C97, we fit the observed data points with an arc-tangent function, which has a 
minimal number of free parameters while adequately describing most galaxy rotation curves 
over the range probed by Ha observations. Specifically, we fit the parameters Vq, tq, rt, and 
K;irc of the relation 

V{ri) = Vo + -Kirc arctan (H^] (1) 

by minimizing x^=J2i [ V {^i)-Vohs{ri) ]'^/<jf, where Vq is the systemic velocity, tq is the spatial 
center of the spectrum, is a turn-over radius where the rotation curve goes from steadily 
rising to fiat, Vdrc is the asymptotic circular velocity, Kbs('^i) is the intensity weighted velocity 
centroid at pixel r^, and cXi is its uncertainty. Note that we do not force vq to correspond to 
the photometric ce nter of the galaxy. We use a Levenberg-Marquardt minimization routine 



(jPress et. al.ll 19921 ). with initial guesses made by eye, to obtain the best-fit parameters and an 
error covariance matrix. Since galaxy disks have non-circular motions at the 10 — 20 kms~^ 
level and we do not want the fit dominated by the high SNR data points at the inner parts of 
the rotation curve, we add lOkms"^ in quadrature to all the velocity centroid uncertainties. 
Changing lOkms"^ to 20 kms~^ makes a negligible difference in the best-fitting parameters. 
Changing 10kms~^ to Okms"^ makes little difference in most but not all cases. The typical 
X^/d.o.f. is 0.1 — 1.0 when lOkms"^ is added to the centroid uncertainty and 1 — 2 without 
this addition. 

Rotation curves are assigned flags indicating how well the data sample the flat part of 
the rotation curve. Flag-1 indicates that there are data along the fiat part of the rotation 
curve. Flag-2 indicates that there are data at the turn-over radius, but not beyond. Flag-3 
indicates that an arc-tangent function adequately describes the data but the rotation curve 
is still rising at the last measured data points. We assign flag-4 to all galaxies with rotation 
curves either cannot be fit by an arc-tangent function, or have no detectable Ha. Although 
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the flags are assigned by eye, the flag correlates well with the difference between the fit 
parameter V^irc and the velocity (Knd) at the last measured point of the rotation curve. Out 
of 234 galaxies, there are 64 flag-4, 50 flag-3, 57 flag-2, and 63 flag-1 rotation curves. 

Figure [6] shows rotation curves and best-fit arc-tangent functions for 16 galaxies with 
flags 1 through 4; the four galaxies whose 2-D spectra appear in Figure H] are shown in the 
right-hand column. Vertical line segments in each panel show the lOth-percentile, median, 
and 90th-percentile values of the la error bars on the observed velocity centroids; typical 
values are 4.5 kms~^, 5.6 kms~^, and 11.0 kms~^. The arc-tangent curves provide good 
descriptions of the overall shape of the rotation curve in all cases, but some rotation curves 
show "bumps and wiggles" associated with non-circular motions. As discussed below, our 
primary measure of galaxy rotation speed is the value Vgo of the arc-tangent fit at the radius 
containing 80% of the z-band flux. This point is marked by an open square in each panel. 
The horizontal lines have a full width of twice 2.2Rd, where Rd is the disk exponential scale 
length (see §4.3). 



4.2. Rotation Speed Definitions 

For TF purposes, we want to characterize the fitted rotation curve by a single rotation 
velocity. While the asymptotic circular speed Vdrc seems the obvious quantity to use when 
the rotation curve is fully constrained, the fitted value can vastly overestimate the true 
rotation speed when the observed data points do not reach the flat portion of the rotation 
curve. The velocity V^nd evaluated from the arc-tangent function at the location of the 
last data point always provides a well constrained rotation speed, with no extrapolation of 
the model fit, and it maximizes use of the data on each individual galaxy rotation curve. 
However, Vend is difficult to model theoretically, because the spatial extent of the Ha data 
varies from galaxy to galaxy. A common choice for optical rotation curves is the rotation 
speed V2.2 at 2.2 disk scale lengths, where the rotation curve of a self-gravitating exponential 
disk would peak. We used this measure in our earlier paper for the disk-dominated subset 
of the galaxy sample (P05). However, for galaxies with significant bulges, the value of Rd is 
sensitive to the degeneracies of bulge-disk decomposition, so we do not want to adopt V2.2 
as our primary velocity measure for the full sample. 

After considering a variety of options, we chose to evaluate the arc-tangent function at 
a radius (-Rso) containing 80% of the i-band flux. This velocity measure, which we refer to as 
Vgo, allows a relatively straightforward comparison to galaxy formation theories. For a pure 
exponential disk, Rso is 3.03Rd, but Rso/Rd is smaller for galaxies with significant bulges. 
The empirical logic of choosing R^q is evident from Figure O most of our rotation curves 
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extend close to or beyond Rso, so substantial extrapolation of the rotation curve is rarely 
required, and it is far enough out to be close to Vdrc in niost cases. 

Figure [7^ makes this point quantitatively, showing the cumulative distribution of 2.2R(i 
and the radii containing 60%, 70%, 80%, and 90% of the i-band flux divided by the radius 
Rend of the outermost Ha data point. Half of the galaxies have data extending beyond i?80- 
The radii 2.2Rd, Rm, and Rjo would under-utilize the data for a large fraction of the sample 
and be further from Vdrc, while evaluating the velocity at Rgo would require extrapolation in 
most cases. Figure [7b shows a generally good correlation between Vgo and Knd, with a hand- 
ful of fiag-3 galaxies having Vso significantly greater than Vend, and some fiag-1,2 galaxies with 
extended rotation curves having Ver^d > Vm- Ou r choic e of Vm as a rotation speed definition 
i s very similar to that of iPersic. Salucci. &: Stell (Il996l ) and ICatinella. Giovanelli. &: Haynes 
who evaluate the rotat i on cu rve at the radius containing 83% of total flux (see 
Catinella. Haynes. fc Giovanelli! (120051 ) for further discussion). In Appendix A, we present 
TF fits for the alternative velocity definitions V2.2 and Vend- 



4.3. Bulge-Disk Decomposition 



Bulge and disk parameters of our target galaxies are needed for three reasons: to correct 
the measured velocities (Vgo) for disk inclination, to correct absolute magnitudes for internal 
extinction by dust in the disk, and to obtain structural quantities that can be tested for 
correlations with TF residuals. The DR2 isophotal axis ratio may be affected by the presence 
of a bulge and spiral arms. Therefore, refitting the disk, separate from the bulge, provides 
more accurate disk axis ratios for the inclination correction of the velocities and the internal 
disk extinction correction. The bulge-disk decomposition also allows measurements of the 
disk size and surface brightness, which are possible third parameters in the TF relation. 



We fit a Sersic profile (ISersid Il968l ) to the bulge and an exponential profi le to the 
disk, using the two-dimensional profile fitting program GALFIT ( iPeng et al.ll2002l ). We run 
GALFIT on the g, r, and z-band images of our galaxies, holding the disk profile fixed as 
an exponential disk and leaving the bulge Sersic index ri as a free parameter. The Sersic 
index describes the central concentration of the light profile, with n=l for an exponential 
disk and n=4 for a DeVaucouleurs profile. We run GALFIT setting the initial bulge Sersic 
concentration equal to n=2 and visually inspect the results. The is sensitive to small-scale 
asymmetries and variations in the galaxy profile that are not two-dimensionally symmetric 
(i.e. HII regions, spiral arms, and bars), so visual inspection of the fit results is required 
to ensure that GALFIT does not drift to a local minimum in an un-realistic location 
of parameter space. Masking out strongly asymmetric features and re-running GALFIT 
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resolves these problems when they arise. 

Figure [S] shows examples of 2-D symmetrical Sersic profile fits to our flag-1,2,3 galaxies. 
The top and middle rows show typical results for galaxies that do not have prominent spiral 
arms. The bottom panel shows a result for a galaxy that has prominent spiral arms. Since 
our galaxies are selected to be inclined, we rarely find prominent bars in the images. 

Figure M shows 1-D profiles along the major axis of the data, model, and residual 
images for galaxies shown in Figure [HI From Figure M "we can conclude that our best-fit 
models accurately describe the galaxy surface brightness profiles down to ~25 mag/arcsec^ 
or fainter, except for small scale variations. The bulge-to-disk flux ratio, disk axis ratio, 
and disk position angle are generally robust outcomes of the fitting procedure, converging 
quickly with little sensitivity to the initial parameter choice. The bulge radius, Sersic index, 
and position angle are subject to more uncertainty, in part because the bulge component 
is often similar in size to the seeing FWHM. Figure [10] shows that approximately 20% of 
our galaxies with usable rotation curves have significant bulges, where early-type galaxies 
have concentrations around 2.6. Figure [TU] also shows that galaxies with lower concentrations 
tend to have higher disk-to-total flux fractions. A visual inspection of the fiag-1,2,3,4 galaxies 
confirms this trend, where the fiag-3,4 galaxies tend to have early morphological types. The 
fiag-1 galaxies are predominantly Sb and Sc Hubble types, whereas the fiag-4 galaxies are 
predominantly Sa with a few Sc galaxies. The flag-2,3 galaxies are mixed between Sa, Sb, 
and Sc types. 



4.4. Internal Extinction Corrections 



Dust in galaxy disks absorbs a larger fraction of the disk light in edge-on directions. 
Therefore, it is standard practice to apply internal extinction corrections to lum inosities used 
for the TF relation. We follow this practice here and use the iTuUy et al.l (Il998l ) formulation 
adapted to our bands. We apply internal extinction corrections to the disk, whi le assuming 



that t he bulge has no correction to the face-on value because it has little dust. iTuUy et al. 



( 119981 ) provide prescriptions for the internal extinction corrections in the Johnson B (438 
nm), R (641nm), I (798 nm), and K' filters as a function of the galaxy inclination and the 
absolute magnitude in that band. We convert the SDSS g (469 nm), r (617nm), i (748nm), 
and 2:-ban d (893 nm) a b solute magnitudes to Johnson magn itudes using the co nversions in 
Table 7 of [Smith et al.l (|2002[ ). then linearly interpolate the [TuUy et al.[ (119981 ) corrections 
back from the Johnson central wavelengths to the SDSS central wavelengths as listed above. 
This internal extinction correction is then applied to the disk g, r, i, and z-band disk fluxes, 
which are then added to the un-corrected bulge fluxes. 
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Figure [TT] shows histograms of the total internal extinction correction in each band. 
Typical values of the (7-band internal extinction range from 0.2 to 0.8 magnitudes. The 
median internal extinction correction decreases from 0.6 in the gf-band to 0.2 in the z-band. 
Figure [T2] shows the (7-band TF relation with and without the internal extinction correction. 
Solid lines show fits using the maximum likelihood procedure discussed in §5 below. The slope 
of the extinction-corrected TF relation is similar to the uncorrected slope (—5.5 ± 0.2 versus 
— 5.2±0.2), though of course the extinction correction increases the zero point by the average 
value of the extinction. The lower panels of Figure [T^] show that the internal extinction 
correction successfully removes a weak trend of TF residual with axis ratio, so it appears 
to be a valid and useful correction on average, even if it is uncertain on a galaxy by galaxy 
basis. We assume, somewhat arbitrarily, that the inclination correction uncertainty is 1/3 
of the correction itself. Results in Appendix A show that, with this assumed observational 
uncertainty, the internal extinction corrections do not change the inferred intrinsic scatter 
a. All magnitudes discussed in the remainder of the paper have this correction applied. 



We are ultimately interested in the slope, intercept, and intrinsic scatter of the TF 
relation. To estimate these, we need an accurate characterization of the observational errors 
of our data. This is especially important for estimating intrinsic scatter, since the total 
scatter about the fitted relation is, roughly speaking, the quadrature sum of the intrinsic 
scatter and the observational errors. 

The TF relation is the correlation of the logarithm of the rotation velocity with the 
absolute magnitude. We define the logarithm of the rotation velocity to be 



where i is the disk inclination angle inferred from the bulge-disk decomposition procedure 
described in §4.3, and the (1 + z) term accounts for cosmological broadening. We relate the 
inclination to the observed axis ratio using the equation 



4.5. Error Budget 




(2) 




(3) 



where h/ais the z-band disk axis ratio determined usin g GALFIT, and 0.19 t a kes in to account 
the finite thickness of the disk (see the discussion in iHaynes fc Giovanellil (119841 )): varying 
0.19 over the range 0.10—0.25 causes the inclination correction to vary by < 1.5%. We define 
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the absolute magnitude to be 

Ma = mpetrcA - 5 logio(f/ig/(l + ^)) + 5 - A\ - w + K,,„,^, (4) 

where A denotes the band {g, r, i, or z), mpetro.A is the SDSS Petrosian magnitude in the 
A-band, A\ is the internal extinction correction described in ^4.4, is the correction for 



Milky Way extinction taken from the SPS S database (ba s ed onlSchlegel. Finkbeiner. fc Davis 



19981 ). -ft'corr,A IS the K-correction from iBlanton et al.l (l2003bl ). and dig is the luminosity 



distance to the galaxy from the Local Group. As noted earlier, we compute distances using 
the SDSS heliocentric redshifts c orrected to the rest frame of the Local Group barycenter 



( lYahil. Tammann. Sandagelll977l ). assuming a cosmological model with fi^ = 0.3, ilx = 0.7, 
and h = 0.7. There are six galaxies (J233152.99-004934.4, J144418.37+000238.5, J095555.07- 
001125.0, J112346.06-010559.4, J203523.80-061437.9, J001006.62-002609.6) that have poorly 
estimated Petrosian magnitudes in the DR2 photometric pipeline. For these galaxies, we use 
the total Sersic magnitude measured using GALFIT. 

The uncertainty in the absolute magnitude is calculated using standard propagation of 
errors, 



{6Mx)^ = (5mpetro,A)' + (5 X 0.434^^ + (^5 x 0.434^^ + {6A 



i \2 

\) 5 



(5) 



where 0.434 converts from natural to base-10 logarithms. Uncertainties in apparent mag- 
nitude 5mpetro,A and redshift 6z are taken from the DR2 database, though we impose a 
minimum redshift error of cSz = 30 kms^^. The uncertainty in the Milky Way foreground 
extinction, A^^ , is assumed to be negligible. We incorporate a distance uncertainty corre- 
sponding to SVpec = 300 kms~^ when calculating disk scale length and lu minosity uncertain- 



ties, to account for the typical amplitude of small scale peculiar velocities (IStrauss fc Willick 



19951 ). We assume the uncertainty due to the internal extinction correction is one third of 
the calculated value, 6A\ = A\/3, and we ignore any (much smaller) uncertainty due to the 
uncertainty in the GALFIT determined i-band axis ratio. In practice, peculiar velocity and 
internal extinction uncertainties completely dominate the uncertainty in Mx- The 10-th, me- 
dian, and 90-th percentile values of the peculiar velocity uncertainty, 5 x 0.434^^^, are 0.049, 
0.080, and 0.118, respectively. The 10-th, median, and 90-th percentile total i-band internal 
extinction correction uncertainties are 0.041 mag, 0.088 mag, and 0.16 mag, respectively. 

The velocity width has uncertainties due to the inclination correction, measurement of 
Vso from the rotation curve, and possible systematic uncertainties due to slit misalignment. 
As described in §2, we removed six galaxies from our sample for which misalignment could 
cause a >10% change in the circular velocity; we do not correct six other galaxies with slight 
slit misalignments for which the correction is <10%. The uncertainty in the inclination 
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corrected velocity width (Vsci) due to the inchnation uncertainty and uncertainties in the 
measurement of Vgo is 

SVsoA' _ (SVsoy ( 5{h/a){h/a) V 

^80. ; Wso ; V 1 - ip/af ) ■ ^ ^ 

The first term dominates for essentially all of the galaxies. Using the covariance matrix 
returned by the Lcvenberg-Marquardt method, the uncertainty of the arc-tangent function 
measured at i?80 is related to the uncertainty in the arc-tangent function parameters via the 
equation 

'•y«,-[i_^c., j. (7) 

The indices run from 1 to 4 for the four arc-tangent function parameters and aj. The 
axis ratio is determined from the disk component after running GALFIT on the i-band 
corrected frames, with the uncertainty reported by GALFIT. The lOth-percentile, median, 
and 90th-percentile la uncertainties in log Vso,i are 0.006, 0.011, and 0.044 logio(km/sec) 
respectively. For a TF slope of —6.0 mag / logio (km/sec), these uncertainties correspond 
to 0.034, 0.063, and 0.264 mag respectively. A typical uncertainty for our program galaxies 
is thus ~ 0.063 mag from the Vso.i uncertainty, ~ 0.088 mag from the internal extinction 
correction uncertainty, and ~ 0.080 mag from the peculiar velocity, summing in quadrature 
to 0.13 mag. This is much smaller than the intrinsic scatter measured in §5.2. However, the 
variation from galaxy to galaxy is large. For the rest of the paper we will use inclination 
corrected velocities and drop the subscript i from Vso,i- Table 2 summarizes the photometric 
parameters and their errors, and Table 3 summarizes the velocity- width measurements and 
their errors. 



5. The Tully-Fisher Relation 

5.1. Modeling the TF Relation 

We use a maximum likelihood method to estimate the slope a, intercept b, and intrinsic 
scatter a of the TF relation. For the "forward" relation, the independent and dependent 
variables are x — r] = log( sin](i+z) ) y ~ respectively. We assume that the intrin- 
sic scatter is Gaussian in form, so that the probability that galaxy i has a true absolute 
magnitude yi given a true velocity width Xi is 

P{yi\xi) = (27rcr^)~^/2exp 



-{yi 



2(72 



(8) 
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where 

yl = a{xi -Xo) + b (9) 

is the mean value expected for a linear TF relation with parameters a and b. The value of 
Xq is chosen so that there is no correlation (or very little correlation) between the statistical 
errors in a and b. For the "inverse" TF relation, we adopt the same model but with x = 
M\ and y — rj. The forward and inverse relations correspond to different assumptions about 
which physical parameter is "primary" : in the forward relation, the absolute magnitude has 
Gaussian scatter about a mean value determined by the linewidth, and the reverse holds for 
the inverse relation. 

We assume that the observational estimates of f j and yi are Gaussian distributed about 
the true values of Xi, yi, i.e.. 



P{yi\yi) = (27rcTjJ ^/^exp 

and 

p{xi\xi) = (27rcT^ J^^/^exp 



2(7 



x.i 



(10) 



where a^^i-, and ay^i are the measurement uncertainties determined as described in §4.5. The 
values of a, 6, and a can then be estimated by maximizing the log-likelihood 

\a.L ^^\Yip{yi\xi,a^^i,ay^i). (12) 

i 

The likelihood for an individual data point can be written 

/oo 
dyiP{yi\yi)p{yi\xi), 
-oo 



-oo 

fOO 



/oo POO 
dyip{yi\yi) / dxip{yi\xi)p{xi\xi). (13) 
-oo J —oo 

We set ^(xjlfj) = p{xi\x.j), implicitly assuming a fiat prior for j9(a;j) over the (typically narrow) 
range of values allowed by the uncertainty in Xi. Substituting equations (8)-(ll) and (13) 
into (12) and simplifying yields the expression 

lBL = -l^ln(a^-^a^ + aV^)-^ '^■:'f' + ';^ + constant. (14) 

Most previous TF studies have implicitly set a to zero when finding the best-fit values of a 
and 6, then estimated a after the fact from the difference between the observed scatter and the 
estimated contribution of observational errors. In the presence of non-zero intrinsic scatter. 
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this method assigns too much weight to the data points with the smallest observational errors, 
yielding statically non-optimal estimates of a and h and underestimates of their uncertainties. 
For example, two data points with near- zero errors would completely dominate the fit, while 
in fact they shoul d be weighted by 1/cr^. Further discussion of these issues can be found in 



D'Agostinil (120051 ). 



The equation 5hiL/5hih = (see eq. 14) can be solved analytically given values of a, a, 
and the input data {xi^Vi^o-x^hO-y^i). We determine the maximum likelihood parameters by 
performing a grid search in a and a, finding the best-fit value of h (analytically) for each (a, 
a) combination, then choosing the (a, 6, a) combination that maximizes InL. We determine 
the la errors on a, 6, and a by repeating this procedure for 100 bootstrap subsamples of 
the full data set, taking the dispersion among the bootstrap estimates as the uncertainty in 
the parameter. We choose the values of Xq in equation (9) so that there is essentially no 
covariance between the error in a and the error in 6; specifically we choose the value of xq so 
that fixing a to a value of ila from its best-fit value does not change the best-fit value of h. 
The resulting zero-points in the r, i, and z TF fits are, respectively, rjo = 2.22, 2.22, 2.22, 
2.23 for the forward fit and Mxo = -20.607, -21.107, -21.327, and -21.400 for the inverse 
fit. 



5.2. The TF Relation in the SDSS Bands 

Figure [T3k shows the z-band TF relation with 162 flag-1,2,3 data points. The data show 
a clear linear trend of with rj and an intrinsic scatter larger than the error bars. The 
best-fit parameters for the forward relation are a = —6.32 ± 0.22, b = —21.390 ± 0.035 (at 
Vgo = 166.0 kms~^), and intrinsic scatter a = 0.423 ±0.035@ Table H] lists these parameters 
and the parameters Oinv, &inv, cTinv of the inverse fit. When comparing the inverse TF to the 
forward TF relation in the figures and text, we invert the inverse relation and refer to the 
slope and scatter by a' = l/oinv and a' = cTinv/ctinv, and quote the intercept b' as the value 
of Mx that corresponds to the zero-point rjo used in the forward fit. The quantities a', b', a' 
have the same units as a, b, a and can be directly compared. For the flag-1,2,3 sample, the 
inverse TF i-band relation a', b', and a' are -7.69 ± 0.29, -21.39 ± 0.04, and 0.47 ± 0.04. 
For all bands and samples, the inverse relation has a steeper slope than the forward relation 
but similar scatter (see discussion in §5.3 below). Throughout the paper we use the terms 
"steep" and "shallow" in reference to the slope (more negative is steeper) of the forward TF 



^The units are mag/logio (kms ^) for the slope and mag for the intercept and scatter. For brevity, we 
wiU omit these units in the text. 
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relation (Vgo predicting Ma). 

Figure [T^.c shows the z-band TF relation when less spatially extended rotation curves 
(fiag-2,3) are removed. As fiag-3 and then fiag-2 galaxies are removed, the slope becomes 
slightly shallower, changing from —6.32 ±0.22 to —6.02 ±0.32 to —5.94 ±0.44. The change is 
at the ~ la level, and it is driven largely by the change in the sample luminosity distribution, 
as there are fewer fiag-1,2 galaxies at low luminosity. Similar trends are seen in g, r, and 
z. The intrinsic scatter declines by a statistically insignificant amount, from 0.42±0.04 to 
0.39±0.05 from fiag-1,2, 3 to fiag-1 only. The constancy of a implies that our fitting procedure 
gives a reasonable estimate of the uncertainty in Vgo even for rotation curves that are rising at 
the outermost point. Given the insensitivity of the results to including fiag-2,3 galaxies, we 
will use the full sample henceforth, for improved statistics and greater sample completeness. 
We list TF parameters for the fiag-1,2 sample in Table E] (see Appendix A). 

Figure [T3] shows the TF relation for the g, r, i and 2;-bands, with parameters fit to the 
data as outlined in §5.1. The forward TF slope increases towards the redder bands, from 
—5.48 to —6.59, a change much larger than the statistical slope uncertainties (typically 0.2). 
This is the expected trend, as fainter galaxies typically have bluer colors. Trends in the 
intercept are not particularly meaningful, as they depend on the adopted (AB) magnitude 
system and stellar mass-to-light ratios in different wavebands. The inverse TF slopes are 
always steeper than the forward slopes, and they show the same trend with wavelength. 
The intrinsic scatter is slightly larger in the ^^-band, and roughly constant in r, i, and 
z. This constancy implies that the extinction corrections are not an important source of 
scatter. It further suggests that the intrinsic scatter is not dominated by variations in stellar 
populations, a point we will return to in §6 below. Excluding fiag-3 and fiag-2 galaxies in 
the g, r, and 2;-bands shows the same trends seen in Figure [131 TF fits using velocity width 
definitions V2.2 and Vend are discussed in Appendix A. 

As shown in §2 (cf. Figured]), the Mj. distribution of our sample is approximately flat 
in the range —18.5 > Mr > —22, but not exactly so. In order to provide a well defined 
target for galaxy formation theories, we have also computed TF parameters after weighting 
galaxies to obtain the results expected for a truly flat Mr distribution. Using the plotted 
distribution of galaxies in 0.5-magnitude bins, we give each galaxy a weight Wi defined as 
the ratio of the number of galaxies in the most populated bin to the number of galaxies in 
its bin. These weights are incorporated into the fitting procedure by multiplying each term 
inside the sums of equation (14) by Wi, as though we had observed the galaxy Wi times {wi 
goes inside the In of the first term). The weights vary from 1 to 5, depending on the bin. 
We do not include the few galaxies with Mr > —18.5 and Mr < —22. The results are quoted 
in Table [61 Overall, the parameters change very little when weighted by our observed Mr 
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distribution. We have carried out a similar weighting experiment to examine the effect of 
exactly reproducing the g-r color distribution of the DR4 "control" sample shown by the 
dotted curve in Figure [3^. The change in the TF parameters (not listed) is much smaller 
than our statistical errors. 

Our sample has absolute magnitude cuts that can create Malmquist-type biases in TF 
estimates. The explicit absolute magnitude-redshift cuts outlined in §2, and the flag-1,2,3 
absolute magnitude distribution seen in Figure [T], may cause a shallowing of a TF slope and 
reduction of the intrinsic scatter when these cuts are applied to data having non-zero scatter. 
Random errors in galaxy distances and scatter across selection boundaries, due to pecuhar 
velocities typically of the order 300 kms~^, can also induce Malmquist-type biases in TF 
estimates. Appendix B describes a Monte Carlo experiment aimed at testing for such biases 
inherent in our sample, assuming a power-law TF relation with Gaussian intrinsic scatter. 
In short, generating a Monte Carlo sample with an input a = -7.35, b = -21.27, and a = 
0.58, then applying the selection criteria outlined above, yields a best-fit forward relation 
of a = -6.32, b = -21.39, and a = 0.42 which is close to the observed i-band TF relation 
(Table 4). The difference between the input and best-fit forward relations is several times 
our statistical uncertainties. The difference between the inverse relations is small. In this 
paper we present a TF relation for a sample of galaxies with well defined selection criteria, 
enabling a theory to reproduce our sample, using our selection criteria, in an attempt to 
model our measured TF relation. Table 4 represents the TF parameters measured with our 
selection criteria. 



5.3. Intrinsic Scatter 

One of our principal objectives is a robust estimate of the intrinsic scatter of the TF 
relation for a broadly selected sample of galaxies. Previous observatio nal studies have show n 



a very wide range of estimates, with some as low as 0.1 to 0.15 mags (iBernstein et al.lll994l ). 
For a large sample of late-type spirals, with similar analysis methods to those used here, 
C97 reports 0.46 mag of total (intrinsic + observational) r-band scatter. K02 calculates 
and - R-band intrinsic scatter of .4 mag for the broadly selected Nearby Field Galaxy Sur- 



vey (jJansen fc Kannapparu l200ll ) and Ursa Major samples, after the selection criteria are 
matched. However, VOl reports a much smaller scatter, 0.15 — 0.18 mag in i?-band, for 
galaxies in the Ursa Major cluster, using the fiat part of the HI rotation curve as the circular 
velocity measure (see Table 6 of VOl). This estimate arises after subtracting 0.17 mag of 
scatter attributed to the depth of the cluster, and changes to the sample or velocity width 
definition can boost the i?-band scatter in the VOl analysis as high as 0.54 mag. Observa- 
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tional esti mates from other studie s span the full range of values quoted above (see e.g., the 



review by IStrauss &: Willickl Il995l ) 



Many papers do not clearly distinguish the intrinsic scatter from the total scatter. In 
part, the observational errors themselves can be difficult to estimate for nearby samples in 
which distance errors dominate and are difficult to define precisely. In addition, previous 
studies have generally not estimated the intrinsic scatter as part of the TF fitting procedure 
as we have done, but have rather subtracted a "typical" observational error in quadrature 
from the total scatter. 

A crucial feature of our sample is that the typical observational errors are smaller than 
the intrinsic scatter because we are measuring galaxies with V > 5000 kms~^ and have 
photometry that allows precise estimates of magnitudes and axis ratios. The most debatable 
elements of our observational error budget (see §4.5) are the 300 kms~^ peculiar velocity error 
and the internal extinction error of 1/3 the internal extinction correction. These typically 
contribute 0.080 and 0.088 mag to the absolute magnitude error. If we increase the assumed 
peculiar velocity uncertainty from 300 to 500 kms~^, per galaxy, our estimate of the i-band 
intrinsic scatter decreases from 0.42 to 0.40 mag. If we set the assumed internal extinction 
error to zero, the estimated intrinsic scatter increases from 0.42 to 0.43 mag. Thus, our 
estimate of a is relatively insensitive to the uncertain elements of our error budget. Even if 
we take the drastic step of setting all of our observational errors to zero, the best-fit intrinsic 
scatter only increases to 0.48 mag. 

When computing inclination corrected velocities, we take the uncertainty in sini from 
the uncertainty in the axis ratio of the GALFIT bulge-disk decomposition. We thus implicitly 
assume that the disk is adequately described by a circularly symmetric exponential, and disk 
ellipticities by d efiniti on contribute to the intrinsic rather than the observational scatter. 



Zaritsky fc Rixl (jl995l ). from a study of face-on galaxies, estimate typical ellipticities of ~ 
0.05 for the gravitational potential in the disk plane, which would cause ~ 0.15 mag of TF 
scatter from a combination of inclination correction errors and non-circular motions (see also 



Franx fc de Zeeuwlll992l : iRvdenI 120051 ) 



Since we have broader sample selection than most previous studies, it is interesting to 
ask whether outliers have an important impact on the TF scatter (or other TF parameters). 
Figure [T5h shows the i-band TF relation with the seven top contributors to marked. From 
largest to smallest Ax^ they are J021941. 13-001520.4, J235106.25+010324.0, J124428.85- 
002710.5, J204913.40+001931.0, J005650.61+002047.1, J013142.14-005559.9, and J235607.82+003258.1. 
None of these galaxies show obvious rotation curve anomalies, but some are morphologically 
unusual. J021941. 13-001520.4 has extended and asymmetrical spiral arms, and its GAL- 
FIT axis ratio is 0.51, giving it a substantial inclination correction. It is therefore possible 
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that the low rotation speed of this galaxy arises from an inaccurate inclination correction. 
The second galaxy, J235106. 25+010324.0, has normal morphology, as does the fifth galaxy 
J005650.61+002047.1. The third galaxy, J124428.85-002710.5, has a prominent central point 
source that produces a persistent residual in the bulge-disk decomposition. However, it is 
under-luminous relative to the mean TF relation, so an AGN contribution cannot explain the 
anomaly. The fourth galaxy, J204913. 40-1-001931.0, has a prominent dust lane, as does the 
sixth, J013142. 14-005559.9. Finally J235607.82+003258.1 has a large ring ~ 28 kpc in size 
centered on the galaxy. However, there are other galaxies with dust lanes or morphological 
asymmetries that are not TF outliers, so we do not think that the intrinsic scatter of the 
sample is simply driven by rare "oddballs". 

Figure \T5b shows the impact of removing, in succession, the data points with the lowest 
likelihood (beginning with the seven shown in Figure [T^). then refitting the TF relation. 
Filled circles show the estimated forward scatter as data points are removed, while open 
circles show the corresponding quantity cr' = (Xinv/flinv of the inverse relation. The inset 
panel shows the best-fit slopes a and a' = l/ainv Removing the first data point, J021941.13- 
001520.4, produces a noticeable drop in a, from 0.42 to 0.40 mag. Since our estimated 
inclination for this galaxy could be inaccurate for the reasons mentioned above, we think 
there is a reasonable case for lowering our estimated values of the intrinsic TF scatter by 
~ 0.02 mag (less than our statistical uncertainty of ~ 0.035 mag). However, removing sub- 
sequent data points produces only a steady, approximately linear decrease of the estimated 
intrinsic scatter. If the TF residuals were drawn from a Gaussian of width 0.4 mag, then the 
estimated scatter should show an approximately linear decrease to zero at Nremoved = 160 (a 
claim we have tested with Monte Carlo experiments). The fact that our estimated scatter 
goes to zero at Nremovcd ~ 90 suggests that we may have overestimated the observational 
errors. If we set the internal extinction correction and peculiar velocity uncertainties to zero 
and repeat the above procedure, the estimated scatter of the forward TF relation approaches 

zero at Nremovcd ~ 130. 

Closely related to the role of outliers is the question of the residual distribution. Fig- 
ure [16] plots the histograms of Ay /ax, where Ay is the deviation of each galaxy from the 
mean forward (solid) or inverse (dotted) relation and 0"^= [c"^ + Cy+o^cr^]^^^ is the quadrature 
sum of the intrinsic scatter and the galaxy's observational error. With the exception of the 
one largest outlier in the forward and inverse relations, both histograms are approximately 
Gaussian in form. (The largest outlier galaxy has Ay far = —4.2.) This is further evidence 
that the intrinsic scatter of our sample is not driven by rare outliers. 

Although the differing slopes of the forward and inverse relation (specifically, the fact 
that a 7^ l/ainv) makes them appear superficially different. Figure [16] shows that they both 
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lead to nearly Gaussian residual distributions. The forward and inverse relations with Gaus- 
sian intrinsic scatter are equally good, and essentially equivalent, descriptions of the two- 
dimensional distribution of galaxies in the {M\,rj) plane, over the range covered by our 
sample. Theoretical models of the galaxy population should explain this full distribution, 
not just a slope and intercept whose values necessarily depend on the fitting procedure. 
Similarly, there is no particular virtue to "orthogonal" fitting procedures that treat the two 
observables symmetrically - they would have intermediate slopes, but they would presum- 
ably lead to a similar two-dimensional distribution once the scatter was properly accounted 
for. 

Our estimate of a = 0.42 mag (for the flag-1,2,3 i-band forward relation) is larger than 
at least some previous estimates of the TF intrinsic scatter. It is interesting to explore 
whether this is a consequence of our broader morphological selection. Figure [T7] compares 
the i-band TF relation of our full sample to a disk-dominated subset that have, according to 
the GALFIT bulge-disk decomposition, i-band disk-to-total flux fractions D/T > 0.9 (see 
P05 for an extensive discussion of this subset). The intrinsic scatter drops io a = 0.36 mag, 
a statistically significant but modest decrease. Selecting disk-dominated systems does not 
yield the small values of the intrinsic scatter found in some previous stu dies. We do not 



see a ny clear evidence for a luminosity dependence of the intrinsic scatter (iGiovanelli et al. 



19971 ). but our sample is relatively small for detecting such an effect. 



Figure [T7b shows the (7-band TF relation of the full sample, with barred, asymmetric, 
and "interacting" galaxies marked separately. Barred galaxies are identified by visual inspec- 
tion of the i-band images, since bars are usually red. We find 17 galaxies with discernible 
bars, but this is an under-estimate of the true barred fraction because bars are difficult to 
detect in highly inclined systems. The asymmetric galaxies are selected by visual inspection 
of the (yf-band images, because warps and subtle interactions cause a boost in star formation 
and bluer colors. We find 27 galaxies with asymmetries. The two "interacting" galaxies 
have close companions. Figure [T7b shows no clear evidence for any of these three classes 
of galaxies to be systematically offset from the full TF relation slope, intercept, or scatter. 
Separate TF fits to the barred and asymmetric subsets yield intrinsic scatter estimates of a 
= 0.29 ± 0.09 and 0.40 ± 0.08, respectively. We conclude that the slightly larger scatter of 
the full sample is not driven by these morphologically distinct populations. 



5.4. Comparison to Previous Studies 



Figure [18] compares our TF results to those of two previous studies: C97, whose analysis 
techniques are similar to our own, and VOl, who presents a comprehensive investigation of 
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spirals in the Ursa Major cluster. C97's sample, shown in Figure [T8b . consists of 304 Sb-Sc 
galaxies selected from the UGC For his data points, we use his preferred definition of rotation 
speed as the amplitude of the arc-tangent function fit to the Ha rotation curve at 2.2 disk 
scale lengths, while for our data points we show VgQ. We ignore slight differences between the 
Gunn r-band and the SDSS r-band. Relative to our sample, C97's is much more strongly 
weighted towards luminous galaxies. Despite these differences, the TF relations agree well 
(solid and dashed lines), with our slope a = —5.96 ± 0.20 slightly shallower than C97's a = 
—6.36 ± 0.22 (C97, Table 4). (If we fit the C97 data with our routines, we obtain a similar 
result.) 

Figure [T8b shows the i?-band TF relation for 29 Ursa Major cluster galaxies from 
vol, with circular velocities defined from the flat portion of rotation curves measured by HI 
synthesis imaging, compared to the r-band TF relation measured here. The VOl TF relation 
is clearly steeper; most of the difference is attributable to the low luminosity galaxies, which 
rotate faster at fixed Mr, though the clump of galaxies at ^ —21.5 is also rotating slightly 
slower. The obvious potential culprit is the velocity width definition, with HI rotation curves 
yielding circular velocities that are systematically higher than Vso for low luminosity galaxies. 
However, using the trends found for the subset of C97 galaxies with both Hct and HI velocity 
widths, P05 conclude that this difference can only account for about half of the difference 
in the TF slopes. Since the number of low luminosity galaxies in VOl's sample is small, the 
impact of this difference in velocity definition could perhaps be enhanced by small number 
statistics. Alternatively, the finite depth of the Ursa Major cluster could contribute to a 
steeper slope if the fainter galaxies happen to lie preferentially on the far side, but this is 
un-likely. The bias-corrected r-band slope, discussed in Appendix B, steepens to a value 
(—7.02) much closer to the VOl -R-band TF sample. Therefore, one can conclude that the 
major differences are a combination of the velocity width definition and sample selection in 
this study. A larger sample of systems with both Ha and HI synthesis rotation curves would 
help shed light on the origin of this difference. 

In Table [7] we present the slopes and intercepts for several TF relations commonly 
used to compare theoretical predictions to observations. The intercepts are not hsted due 
to variations in the assumed value of Hq, and different internal extinction corrections. The 
difference between slopes, in a given band, is several times the typical statistical uncertainties, 
with no obvious trends. As shown above, different velocity width measurements can account 
for some of those differences. K02 is also able to measure slopes for the MAT and CF 
field galaxy samples that are similar to the K02 slopes after applying a consistent internal 
extinction correction, velocity width definition, and fitting method. The intrinsic scatter is 
typically smaller in the studies using HI velocity widths. Those samples are primarily used as 
distance indicators to study peculiar velocity fields, and represent more pruned samples. For 
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example, the small scatter measured by iPierce fc Tullyl (Il992l ) was measured for a pruned 
sample of nearby spiral galaxies having well measured distances. In samples with optical 
velocity widths and careful considerations of biases, the intrinsic scatter appears to be 0.3- 
0.5 magnitudes. Given our sample's broad selection criteria in terms of galaxy Hubble type, 
well defined magnitude and redshift cuts, and further distances we consider our TF relation 
to suffer less peculiar velocity uncertainties and Hubble type selection effects. However, the 
large differences in the slopes may be due to the Malquist-type biases discussed in Appendix 
B. 



6. Residual Correlations 



To understand the physical sources of scatter in the TF relation, we want to investigate 
the correlation of TF residuals with other galaxy properties. For examp le, in the case of 
ellipt icals, the correlation of residuals from the Faber- Jackson relation (IFaber &: Jackson 
19761 ) with galaxy size led to the recognition that ellipticals occupy a fundamental plane 
that l argely corresponds t o the virial relation for the stellar component (iDjorgovski fc Davis 
19871 : iDressler et al.lll987l ). The lack of a similar correlation for disks implies that disk gravity 
does not dominate the ro tation speed at radii used for TF investigations (ICourteau &: Rix 
19991 : iPizagno et al.ll2005l ). If the rotation velocity of disk galaxies is more fundamentally 
related to the stellar mass than to the luminosity, there should be a correlation of TF residual 
with color, which tracks the stellar mass-to-light ratio (K02). F or a theoretical di s cussio n 
of some of these p oints, see IConti. Ryden. fc Weinberg! (120011 ). IShen. Mo. fc Shul (120021 ). 
Dutton et all J2007h . and ICnedin et all J2006h ~ 



The top panels of Figure [11] show the forward and inverse ^f-band TF relations. The 
middle panels show the correlation of the extinction corrected g-r color with rj and Ma. 
The lines show the maximum likelihood best-fit relations, with parameters and bootstrap 
uncertainties listed in Table [51 Point types encode the residual color relative to this mean 
relation, with filled circles, open circles, and triangles showing the reddest, intermediate, and 
bluest 1/3 of the galaxies. The same point type is used for each galaxy in the upper panels, 
and one can see that red galaxies tend to be slightly underluminous in the forward relation. 
The top right panel of Figure [T9l clearly shows that red galaxies tend to rotate faster at fixed 
Mg. The bottom panels plot the residual from the TF relation against the residual from 
the color-?7 or color- Mg relations; solid lines show the maximum likelihood fit to the mean 
correlation of residuals, and dotted lines show bootstrap uncertainties. The residuals are 
correlated, again more clearly for the inverse relation, but there is substantial scatter that 
is large compared to the observational errors. Note that while the bootstrap errors on the 
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best-fit linear slopes are small, the data are not well described by any linear relation with 
zero intrinsic scatter, and the derived slopes therefore depend significantly on the fitting 
procedure. Figure [20] shows similar results for the i-band TF relation. The correlations 
between color and z-band TF residuals are much weaker, being essentially absent in the 
forward relation. 



In the population synthesis modeling of iBell et al.l (120031 ). the mass-to-light ratio of a 



stellar population changes with g-r color roughly as M^jLg (x {g — r)^'^^ in the ^f-band, 
and M^,/ Li (x {g — r)°-^^ in the i-band. For pure self-gravitating disks, the circular velocity 
should correlate with the stellar mass as oc M^, at fixed scale length. Variations of the 
stellar mass-to-light ratio would therefore produce inverse TF residual correlations of the 
form A{g — r) ~ 2Ar](Mg) /1.52 and A{g — r) ^ 2 Ar] (Mi)/ 0.86, and forward TF residuals of 
the form A{g-r) ~ -0AAMg{r])/1.52 and A{g-r) ~ -0.4AMj(r/)/0.86, in the absence of 
other effects. Dashed lines in Figure [19] and Figure [20] show the correlation slopes predicted 
by this simplistic model. 

The residual correlations in the bottom panels of Figure [TH] and Figure [20] have the 
correct sign expected for variations of M*/L with stellar populations. Therefore, we concur 
with the conclusion of K02 based on the NFGS, that these variations account for some 
of the scatter in the TF relation. However, it is clear from the scatter about the mean 
residual correlation, especially in the i-band, that these variations do not account for much 
of the intrinsic scatter. We have tried the experiment of changing rj for each galaxy by an 
amount —Arj predicted from its g-r color using the best-fit slope to the inverse TF residual 
correlations from Figure [19] and Figure [201 then refitting the TF relation. This procedure 
reduces the estimated intrinsic scatter of the inverse relation from a' = 0.073 to 0.057 in the 
^f-band and from a' = 0.061 to 0.057 in the i-band, drops of 22% and 7%, respectively. 

We can also explore the correlations of TF residuals with structural parameters. Fig- 
ure ?? presents the correlation of TF residuals with z-band half-light radius Ri, determined 
from the GALFIT model fits to the z-band images, in the same format as Figure [19] and 
Figure [20] In the forward relation, there is a slight tendency for larger disks at fixed rj to be 
slightly more luminous (upper left panel). This trend leads to a weak correlation between 
the TF residual and the residual from the mean Ri-r] relation, though the scatter is large 
compared to the mean correlation (and to the observational errors). The inverse fits reveal 
no trend of residual Vso with residual Ri at fixed Mj. Dashed lines show the predictions of 
a pure self-gravitating disk model, with V"^ (x 1 / Ri at fixed Mj. 



These weak or absent residual-radius correlations confirm the results of Courteau fc Rix 



([1999]), now with a more broadly selected sample and smaller observational errors per galaxy. 



We found similar results for a disk-dominated subset in P05, using V2.2 as a rotation measure. 
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For the full sample, we have also investigated using the velocity at radii c ontaining 60% to 
90% o f the total i-band flux, again with sirnilar r esult s. As discussed bv ICourteau &: Rix 
(119991 ) ■ and in greater detail by lDutton et al.l (120071 ) and lGnedin et al.l (120061 ) . the absence of 
strong radius residuals imposes strong constraints on the contribution of disk gravity to the 
rotation speed; disk gravity should cause more compact galaxies t o rotate faster, and the im- 
pact of the disk on the inner halo profile (IBlumenthal et al.lll986l : iGnedin et al.ll2004j ) should 
amplify this effect. Explaining the observed lack of correlation requires "sub-maximal" disks, 
and it is more easily accomplished if disks do not produce adia batic contraction of halos (see 
Dutton et allboO?! : ICourteau fc Rb3ll999l : ICnedin et alll2006h . 



Figure 122] examines the dependence of the TF residuals on morphology, as quantified by 
the concentration index Cj = Rgo/R^o, where Rgo and R^q are the radii enclosing 90% and 
50% of the 2-band Petrosian flux, as determined by the SDSS photometric pipelines . An 
index of c ~ 2.6 is often ad opted as the separati on between early-type (high concentration) 
and late-type galaxies (e.g.. IStrateva et al.ll200ll ). We plot early-type and late-type galaxies 
as filled and open circles in the top panels of Figure [22l using the q = 2.6 division, and the 
bottom panels plot q against the TF residuals. On average, early-type galaxies tend to be 
fainter at fixed r] or rotating faster at fixed Mj, though the mean offset is small compared 
to the scatter. The mean residuals for early-type galaxies are < AMj > = 0.14±0.09 and 
< Ar]> = 0.026±0.015, rising to 0.19±0.08 and 0.034±0.012 if we exclude the single strong 
outlier (J021941. 13-001520. 4) discussed in §5.3. This trend agrees with the finding by K02 
that Sa galaxie s are f ainter than later type spirals at fixed rotation speed ( see also G97 and 
Masters et. al.l (120061 )). The sign of the observed trend agrees with expectations from older, 
redder stellar populations of early-type galaxies. However the separation in Figure [221 while 
modest, appears somewhat clearer than the trend with color seen in Figure [20l suggesting a 
direct correlation of rotation speed with morphological structure. Using the GALFIT D/T 
ratio, instead of Cj, gives similar results. 



7. Conclusion 

We have measured the TF relation for a sample of 162 galaxies selected from the SDSS, 
using follow-up spectroscopy to obtain Ha rotation curves. We targeted 234 galaxies from 
the SDSS spectroscopic galaxy catalog that have a roughly fiat absolute magnitude distribu- 
tion in the range —18.5 > > —22. For the purpose of testing galaxy formation models, 
our sample has several advantages relative to most previous TF studies. Our target selection 
is blind to morphology, except for an z-band axis ratio cut of 0.6. Our completeness is high, 
and the galaxies with usable Ha rotation curves have distributions of color and concentration 



-29- 



similar to a large control sample of galaxies selected from DR4 with the same Mr distribu- 
tion but no inclination cut or Ha requirement. Due to the uniform SDSS photometry and 
the selection of galaxies at rcdshifts greater than 5000 kms~^, all sources of observational 
error arc small compared to the estimated intrinsic scatter. The uniform multi-wavelength 
SDSS photometry allows us to investigate the TF relations and the correlations between TF 
residuals, color residuals, and other structural parameters in the g, r, i, and 2;-bands with 
minimal photometric uncertainties. 

We adopt Vso as a measure of the rotation speed and use extinction corrected absohitc 
magnitudes. Wc estimate the slope, intercept, and intrinsic scatter of the TF relations simul- 
taneously, using a maximum likelihood procedure that accounts for individual observational 
errors in absolute magnitude and rotation speed. We measure forward TF slopes between 
—5.5 and — 6.6 mag/logio( kms~^), with typical uncertainty of 0.2, and intercepts between 
—20.7 and —21.5 mag at 77 = 2.22, with typical uncertainty of 0.04. The slope becomes 
systematically steeper with wavelength. Inverse fits always yield steeper slopes, but once 
the effects of Gaussian scatter are included, the forward and inverse fits describe essentially 
the same two-dimensional distribution of data points over the range covered by our data. 
Corrections to these slopes due to Malmquist-type biases are discussed in Appendix B, where 
it is shown that our observed TF relation can be reproduced with a steeper slope and larger 
intrinsic scatter given certain assumptions about selection effects. 

The intrinsic scatter appears to be nearly independent of wavelength or fitting procedure, 
typically a = 0.42 — 0.46 mag, with a higher value (0.54 mag) for the inverse fit in g-hand. 
The distribution of residuals is approximately Gaussian for both the forward and inverse 
relations, and there is no indication of rare outliers inflating the TF scatter, with the possible 
exception of one galaxy whose axis ratio may be incorrectly measured because of spiral arms. 
Omitting this one system reduces the estimated intrinsic scatter by 0.02 mag. The intrinsic 
scatter is slightly smaller for a disk-dominated subset of galaxies, decreasing in the i-band 
from 0.42 mag to 0.36 mag for galaxies having disk-to-total flux ratios greater than 0.9. 
Morphologically asymmetric, barred, or possibly interacting galaxies show no clear evidence 
for offsets from the mean TF relation or for larger scatter, though our statistics for these 
subsets are limited. 

In an attempt to understand the origin of scatter in the TF relation, we study correla- 
tions between TF residuals and other galaxy properties: g — r color, half-light radius, and 
i-band concentration index. The g-hand TF residual shows a significant correlation with 
color, in the sense that bluer galaxies are brighter at fixed f] or (more clearly visible in the 
data) that redder galaxies rotate faster at fixed Mg. However, these correlations are much 
weaker in the i-band, and they are much weaker than predicted for pure self-gravitating 
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disks with the stellar mass-to-light ratios predicted from g — r color. The correlation with 
half-light radius is very weak, with a slight tendency for larger galaxies to be more luminous 
at fixed rj. More concentrated (earlier type) galaxies tend to be slightly fainter at fixed 77 or 
to rotate slightly faster at fixed Mj, but the trend is weak compared to the TF scatter. 

Our TF relation is similar to that found for field spirals by C97, who used similar analysis 
methods. The TF relation found by VOl for Ursa Major spirals, using the flat portion of 
HI rotation curves as a velocity measure, is significantly steeper than ours. This difference 
may partly reflect the difference in velocity width definitions, and Malmquist-type biases as 
described in Appendix B. The residual correlati ons in our sample a re similar to those found 
by K02 using the Nearby Field Galaxy Survey of lJansen et al.l (l2000l ). but we find that stellar 
population effects explain only a small fraction of the scatter in the TF relation. 

The nearly constant intrinsic scatter of the TF relation has interesting implications 
for galaxy formation theories. The scatter is only slightly larger in bluer bands and only 
slightly smaller for a morphologically selected disk-dominated sub-sample of galaxies, and 
the correlation of TF residuals with other properties explains only a small fraction of the 
scatter. The lack of correlation with radius implies that disk gravity has small influence on 
galaxy dynamics at the radii used for TF meas urements; this observational result prov ides a 
strong constraint on galaxy formation models (jPutton et al.l 120071 : iGnedin et al.ll2006l ). The 
weak correlation with color, size, and morphology implies that most of the scatter in the TF 
relation arises from genuine variation in the ratio of dark matter to baryonic matter from 
halo to halo. 

Our results suggest many fruitful directions for future observational investigations, each 
of them comparable in magnitude to the one undertaken here. Absorption-line velocity mea- 
surements, which require longer exposures than our Ha emission line measurements, could 
yield rotation velocities for the remaining 30% of our sample, showing whether galaxies 
without Ha emission follow a systematically different TF relation. HI measurements would 
allow detailed comparison between TF relations defined by Ha and HI velocity widths (see 
C97) for our broadly selected galaxy sample, and they would provide gas masses and gas 
fractions as additional paramet ers for residual correlation studies, whi ch would be valuable 



for testing theoretical models (jPutton et al.l 120071 : iGnedin et al.l |20C 



^. Extending the T F 



relation to brighter and fainter luminosities is of great interest (e.g.. iMcGaugh et al.ll2000l ). 
though complications arise from the growing paucity of disks in the first case and the in- 
creasing importance of non-circular motions in the second. Finally, it is clear that detailed 
investigation of TF residuals would benefit from a much larger sample with similarly broad 
selection criteria and similarly small observational errors, so that one could, for example, 
investigate residual correlations with size, morphology, or environment within restricted re- 
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gions of luminosity-color space. Sin ce large samples of galaxies (> 2000) have been observed 
for TF studies in the past (see, e.g.. lCourteau et al.ll2003l : [Masters et. al.ll2006l and references 
therein), it may be possible to construct such a sample largely from archival data, if the data 
quality can be made sufficiently homogeneous and the effects of sample selection can be suf- 
ficiently well understood. Our data already provide a well defined and highly constraining 
target for theoretical models of galaxy formation. Extended observational studies would 
improve our ability to deduce the essential physics of disk galaxy formation empirically. 
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A. TF Parameters for Alternative Parameters and Sample Cuts 

Our standard TF fits, presented in Table HI use all fiag-1,2,3 galaxies, absolute magni- 
tudes corrected for internal extinction, velocities defined by Vso, and equal contribution of 
each galaxy to the likelihood sums of equation f|T^ . Table [6] presents TF fits using a number 
of variations on parameter definitions and sample cuts or weights. First, we list results for 
the full sample with the usual definitions but weights chosen to approximate a truly fiat 
absolute magnitude distribution in the range —18.5 > Mr > —22 (see §5.2). Next, we list 
results for the full sample with the usual weights but no internal extinction correction (see 
Fig. [T^ . Then we list results for two alternative velocity width definitions: the value of the 
arc-tangent functional fit evaluated at 2.2 disk scale lengths as determined from the i-band 
GALFIT model (V2.2), and at the location of the most distant Ha data point (Vend)- Finally, 
we list the TF fits for the standard definitions and fiag-1,2 galaxies only. Although slopes 
and intercepts change somewhat with these choices, none of our conclusions about the TF 
relation would be qualitatively different for any of these alternative choices. 



B. Malmquist-type Biases in Our TF Sample. 

Applying magnitude cuts to TF samples can lead to biased estimates of the TF pa- 
rameters because there are low and high luminosity galaxies, at fixed velocity, scattered 
above or below the magnitude limits. This can be due to both the intrinsic scatter and 
peculiar velocities scattering galaxies across selection boundaries. Most discussions of these 
"Malmquist"-type biases have focused on apparent magnit ude limited samples, or on the 
systematic bias in the derived peculiar veloci ty field (e.g., iLynden-Bell et al.l Il988l : iGould 
I993I : lTeerikorpilll993l : IStrauss fc Willicklll995l ). Since our selection procedure is quite differ- 
ent from those of most previous TF surveys, we have tested for Malmquist-type biases with 
a simple Monte Carlo experiment. 

We first investigate the affect our absolute magnitude-redshift cuts and fiag-1,2,3 abso- 
lute magnitude distribution will have on our measured TF relation. For our Monte Carlo 
experiment, we generate a sample of 10^ galaxies, evenly sampled in logarithmic velocity 
space, assigning random positions in a space between 4000 kms~^ and 16000 kms~^. We 
then assign absolute magnitudes to each galaxy following an assumed TF relation, having a 
Gaussian intrinsic scatter, represented by 



M = aubiv - Vo) + bub + N{0, cr„fe). 



(Bl) 



where a„f,, but, and aub are the slope, intercept, and intrinsic scatter for an un-biased power- 
law TF relation. We use tjq = 2.22 throughout the procedure. We then apply our absolute 
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magnitude- redshift cuts, outlined in §2, and select galaxies according to our observed flag- 
1,2,3 absolute magnitude distribution, with a final Monte Carlo sample that is 10 times 
larger than our flag-1,2,3 sample, or 1620 galaxies. Generating a Monte Carlo sample with 
d-ub = -7.52, bub = -21.23, and aub = 0.45, and applying our cuts, will have a measured 
TF relation with a slope, intercept, and scatter of -6.86, -21.27, 0.33. This is quite a large 
difference considering our typical statistical uncertainties on the measure slope, intercept, 
and scatter being 0.2, 0.04, and 0.04 mag respectively. The effect on the inverse relation is 
much smaller since the cuts are applied to absolute magnitudes. Figure [23] shows this Monte 
Carlo sample with and without these cuts. It is clear from Figure [23] that the selection 
criteria has excluded high and low luminosity galaxies, shallowing the slope and reducing 
the scatter. 

The measurement errors in SDSS Petrosian magnitudes are generally small, and lumi- 
nosity uncertainties are therefore dominated by line-of-sight peculiar velocities, which we 
have assumed in our analysis to be drawn from a Gaussian of dispersion 300kms~^. Next 
we include peculiar velocities in the Monte Carlo sample by modifying each galaxy's red- 
shift with a peculiar velocity drawn from a 300kms~^ Gaussian. We then apply the same 
absolute-magnitude redshift cuts, while matching the flag-1,2,3 absolute magnitude distribu- 
tion in the shaded area of Figure [T] then measure the TF relation for this sample. Generating 
a Monte Carlo sample with aub = -7.52, bub = -21.23, and aub = 0.45 has a measured slope = 
-6.77, intercept = -21.34, and a = 0.33 mag. This is a relatively small, but non-zero, effect. 
We include the 300kms~^ peculiar velocities in the analysis below. 

We now estimate what power-law TF relation and Gaussian intrinsic scatter would 
be required to reproduce the measured forward TF relations in Table 4, given our selection 
criteria. This is done iteratively by varying aub, Kb, and aub, until the fitted parameters, after 
applying our selection criteria and 300kms~^ peculiar velocities, agree with the measure TF 
relations reported in Table 4. The forward i-band TF relation in Table 4 is reproduced with 
un-biased values of aub = -7.35, bub = -21.27, and aub = 0.58 mag. Comparing these numbers 
to the best-fit values for the forward i-band relation in Table 4, shows that the un-biased 
slope is steeper by ~ 5 sigma, the intercept changed by several sigma, and the un-biased 
intrinsic scatter is larger by ~ 4-5 sigma. Similar results are found for the g, r, and 2;-bands. 
The un-biased r-band TF parameters, that reproduce the best-fit values in Table 4, are aub 
= -7.02, bub = -21.05, and aub = 0.56 mag. 

Table 4 presents the best-fit TF relation to galaxies observed with the selection criteria 
outlined in §2. Although the un-biased parameters discussed here are quite different from the 
best-fit TF parameters, the best-fit parameters have fewer assumptions than the un-biased 
parameters and have well defined selection criteria. Namely, the un-biased power-law TF 
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relation assumes that there is no change in slope at high and low luminosities (contrary to 
K02), and that the best-fit TF parameters in Table 4 are a fit to the data and not a Monte 
Carlo experiment. However, G97 show that these assumptions have a small affect on the un- 
biased parameters. Therefore, the recommended procedure when trying to reproduce the TF 
relation in the SDSS bands would be to apply our selection criteria and absolute magnitude 
distributions to a sample (theoretical or observed), then measure the best-fit parameters. 
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Fig. 1 . — Distribution of r-band absolute magnitudes for the full sample targeted for follow- 
up Ha observations (234 galaxies, full histogram) and the subset of the sample with usable 
(flag-1,2,3) Ha rotation curves (170 galaxies, shaded histogram) 
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Fig. 2. — The color magnitude relation for our sample of galaxies plotted on top of contours 
containing 2 5%, 50%, and 75% o f the total luminosity density in the whole low redshift 
DR2 sample (jBlanton et al.ll2003bl ). The data point type indicates the quality of the galaxy 
rotation curve corresponding to flag-1 (circles), flag-2 (squares), flag-3 (asterisks), and flag-4 
(diamonds). 
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Fig. 3. — Cumulative distributions of the (top) g-r model color and (bottom) i-band con- 
centration index (rgo/rso) for our full sample (short-dashed curves), the subset with usable 
rotation curves (long-dashed curves), and "control" samples with the same Mr distribution 
selected from SDSS DR4, with (dotted) and without (solid) an isophotal axis ratio cut b/a 
< 0.6. 
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Fig. 4. — Examples of flat-fielded, wavelength calibrated, and linearized spectra of four 
galaxies with varying quality flags for their Ha rotation curves. Column (a) shows a flag-1 
spectrum, with data on the flat portion of the rotation curve. Column (b) shows a flag-2 
spectrum, with data extending just past the turn-over radius. Column (c) shows a flag- 
3 spectrum, with a rising rotation curve at the last detectable data point. Column (d) 
shows a flag-4 spectrum, with insufficient Ha for a useful rotation curve. The horizontal 
bars indicate 10 A and the vertical bars indicate 20 arc-seconds. Lower panels show the 
spectra after continuum subtraction. The central "curve" is the Ha emission line used for 
the rotation curve measurement; (a)-(c) also show the [NH] emission lines. 
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Fig. 5. — Comparison of rotation curves measured at MDM (filled triangles) and Calar Alto 
(open triangles). Dotted and solid curves show the arc-tangent function fits to the MDM 
and Calar Alto data, respectively. 
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Fig. 6. — Rotation curves of 16 sample galaxies. The data quality flag is marked in each 
panel, and the four galaxies shown in Figure H] appear in the right hand column. Vertical line 
segments in each panel show the lOth-percentile, median, and 90th-percentile values of the 
la error bars on the observed velocity centroids. Solid curves show the best-fit arc-tangent 
functions. Open squares mark the fitted velocity at the radius enclosing 80% of the i-band 
fliix, which is our standard measure of the galaxy rotation speed. Horizontal bars extend 
from —2.2Rd to +2.2Rd, where Rd is the disk exponential scale length of the z-band GALFIT 
model. 
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Fig. 7.— (Top) Cumulative distributions of the radii that contain 60%, 70%, 80%, and 90% 
of the i-band flux divided by the radius of the last Ha data point. The distribution for 2.2i?d 
is also shown. (Bottom) Comparison of our standard velocity measure Vgo, the value of the 
arc-tangent function evaluated at the radius containing 80% of the i-band flux, to Knd) the 
velocity evaluated at the radius of the last Ha data point. 



-47- 

















/ 


- * ^ . '". 






4 







Fig. 8. — Two-dimensional bulge-disk decomposition of the i-band images of J235656.66- 
005912.3, J094949.62+010533.2, J135946.01+010432.8 (top to bottom). From left to right 
the panels are the data, the GALFIT model, and the residual image given by subtracting 
the GALFIT model from the data. 



-48- 




3^ -40 -20 20 40 -40 -20 20 40 




r{arc-seconds) r(arc-seconds) 

Fig. 9. — Major axis surface brightness profiles (averaged over five pixels in each direction) 
of the three galaxies shown in Figure [HI In the left panels, points show the data, solid curves 
show the GALFIT model, and vertical lines mark the half-light diameter. Right hand panels 
plot the residual counts (model— data); for these galaxies, 10 counts/pixel corresponds to 
roughly 23.4 mag/arcsec^ in the z-band. 
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Fig. 10. — The i-band D/T vs. the concentration index (rgo/rso) for the flag — 1,2,3 
galaxies. 
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Fig. 12. — Effect of the internal extinction correction on the g^-band TF relation. Upper 
panels show the data (points) and best-fit TF relations (sohd lines) without (left) and with 
(right) internal extinction corrections. The dotted line shows the fit from the other panel. 
Lower panels plot TF residual against i-band axis ratio, with solid lines showing the best-fit 
trend and dotted lines the bootstrap-estimated uncertainty. We define t] = log Vgo/ ( km s~^), 
where Vgo is the inclination corrected circular velocity at the radius containing 80% of the 
i-band flux. 
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Fig. 13. — The i-band TF relation for the full galaxy sample (top), flag-1,2 galaxies (middle), 
and fiag-1 galaxies only (bottom). Solid and dotted lines show forward and inverse fits, 
respectively. The forward slope, intercept (at rjo = 2.22), and intrinsic scatter are listed as a, 
b, and cr, while a', b', a' refer to the converted inverse fit parameters. Parameter uncertainties 
are listed in Tables H] and [61 Table |H 
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Fig. 14. — TF relations for the full (flag-1,2,3) sample in the g, r, i, and 2;-bands. The format 
is similar to Figure [131 except that we show the 10-th percentile, median, 90-th percentile 
la errors crosses in each panel instead of plotting errors on each data point. All intercepts 
(6, b') are quoted at the rjo values listed in Table 4. 
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Fig. 15. — (Top) TF relation with the seven largest outliers (in terms of Ax^) marked in 
rank order. (Bottom) Estimated intrinsic scatter of the i-band TF relation after omitting 
the A^rcmovcd data points with the largest A^^. The solid circles represent the forward TF 
relation and the open circles the inverse. The slope as a function of removed data points is 
shown in the inset. 



-55- 




Fig. 16. — Distribution of residuals for the i-band TF relation, with Ay ^ y — y{x), y = Mi 
for the forward (solid histogram) relation and y = rj for the inverse (dotted histogram) 
relation, and the total scatter for a galaxy with observational errors a^, Cy defined by = 
a'^ + ay + {aUxY. The smooth solid curve is a Gaussian of unit dispersion. The TF outlier 
J021941. 13-001520.4 for both the forward and inverse relations is at Ay/uT ~ —4. 
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Fig. 17. — TF relations for morphologically identified subsamples. (Top) Filled circles 
represent galaxies with disk-to-total flux ratios D/T > 0.9; open circles represent other 
sample galaxies. Solid and dotted lines show best-flt i-band TF relations for the full and 
disk-dominated samples, respectively. (Bottom) TF relation in g-hand with asymmetrical, 
barred, and "interacting" galaxies (with nearby neighbors) marked by triangles, squares, and 
asterisks, respectively. The solid line is the best-fit relation for the full data set. 
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Fig. 18. — Comparison to results from C97's field sample and VOl's sample from the Ursa 
Major cluster. (Top) Filled and open circles show galaxies from C97 and from this paper, 
respectively; dotted and solid curves show the corresponding forward TF fits. Velocities 

are Vgo for our points and V2.2 for C97. (Bottom) Same as above, but with filled triangles 
representing the VOl data. Absolute magnitudes are SDSS r for our points and Johnson R 
for vol. Velocities for VOl are Vflat, measured from the flat portion of HI synthesis rotation 
curves. 
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Fig. 19. — Correlation of (/-band TF residuals with g — r color residuals, for the forward 
(left) and inverse (right) relations. In the middle panels, solid lines show the best-fit mean 
trend of (extinction corrected) g — r with rj or Mg, with dotted curves showing the bootstrap 
uncertainties. Filled circles, open circles, and triangles show the reddest, intermediate, and 
bluest 1/3 of the galaxies relative to these mean relations; the same point types are used for 
each galaxy in the upper and lower panels. Bottom panels plot residual from the mean color 
relation against TF residual. Solid and dotted lines show the best-fit residual correlation 
and la uncertainties. Dashed lines show the expected slo pe for pure self-gr avitating disks 
with mass-to-light ratio as a function of color according to iBell et al.l (120031 ). Error crosses 
in each panel show median values of the observational errors. 
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Fig. 20. — Correlation of i-band TF residual with residual from the mean relation between 
z-band half-light radius Ri and Mj or rj. Format is similar to Fig. [13 
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Fig. 21. — Correlation of i-band TF residualfrom the mean relation between i-band half-light 
radius i?j and Mj or 77. The Format is otherwise the same as Figure 19. 
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Fig. 22. — Correlation of i-band TF residual with i-band concentration index Cj = rgo/rsQ. 
Lower panels plot Cj vs. TF residual. In all panels, filled and open circles show "early" and 
"late" galaxy types, respectively, with the division between types at Cj — 2.6. 
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Fig. 23. — A Monte Carlo sample without absolute magnitude cuts (small dots) and with 
our absolute magnitude cuts applied (filled circles). The relation used to generate the sample 
is shown as a solid line, the dotted line is the best-fit to the sample after the cuts have been 
applied. 
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Table 1. Telescope and Spectrograph Parameters 



Parameter 

Telescope 
Spectrograph 
Grating 
Blaze 

Spectral Coverage 

Slit Width 

Spatial Scale 

Dispersion 

Read Noise 

Gain 

Detector 

Net Instrument Eff . 



Calar Alto 



3.5 meter 
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1200 lines/mm 
7100 (1st order) A 
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0.56" /pix 
0.55 A/pix 
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0.865 e-/ADU 
SITe 2000 X 800 CCD 
0.2 



MDM 



2.4 meter 
CCDS 
600 lines/mm 
5875 (2nd order) A 
6500-6994 A 
2" X 5' 
0.41" /pix 
0.41 A/pix 
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2.1 e-/ADU 
Loral 1200x800 CCD 
~ 0.15 



telescope-l-spectrograph-l-detector 



Table 2. Photometric Parameters 
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Table 2 — Continued 
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J234125.92-003943.9'' 


162.37(4.40) 


-19 


94 


,0 


06) 


-20 


80 





06) 


-21 


22 





06) 


-21 


53 





06) 




















J234215.64-003745.9 


160.92(4.40) 


-20 


28 


'0 


14) 


-20 


97 





15) 


-21 


09 





10) 


-21 


26 





08) 





75(0 


02) 





77 


3.84(0.13) 





38(0 


01) 


J234328.26-000148.6 


165.10(4.40) 


-20 


79 





27) 


-21 


35 





22) 


-21 


57 





17) 


-21 


72 





13) 





63(0 


04) 





95 


5.87(0.24) 





25(0 


01) 


J234504.86-001615.1 


103.61(4.36) 


-20 


32 





29) 


-20 


62 





23) 


-20 


71 





19) 


-20 


72 





16) 





37(0 


06) 





98 


7.19(0.33) 





23(0 


01) 


J235106.25+010324.0 


120.59(4.37) 


-22 


20 





28) 


-22 


61 





23) 


-22 


75 





19) 


-22 


44 





14) 





61(0 


06) 





97 


13.95(0.61) 





38(0 


01) 


J235410.09+002258.3^ 


117.04(4.37) 


-21 


47 


,0 


09) 


-22 


35 





09) 


-22 


80 





09) 


-23 


07 





09) 




















J235603.89-000958.6 


116.31(4.37) 


-20 


22 


'0 


20) 


-20 


47 





17) 


-20 


51 





14) 


-20 


45 





12) 





25(0 


04) 





97 


5.57(0.21) 





34(0 


01) 


J235607.82+003258.1 


99.38(4.36) 


-20 


66 





16) 


-21 


29 





14) 


-21 


54 





10) 


-21 


69 





10) 





75(0 


04) 





83 


4.12(0.20) 





70(n 


01) 


J235624.68-001739.6 


110.86(4.37) 


-18 


99 





16) 


-19 


39 





14) 


-19 


51 





12) 


-19 


58 





11) 





39(0 


03) 


1 


00 


3.32(0.15) 





36(0 


01) 


J235634.52+003300.l'^ 


103.38(4.36) 


-17 


96 





09) 


-18 


38 





09) 


-18 


63 





10) 


-18 


82 





12) 




















J235656.66-005912.3^ 


146.95(4.39) 


-19 


78 


,0 


07) 


-20 


57 





07) 


-20 


95 





07) 


-21 


20 





07) 




















J235743.74+003918.6 


103.99(4.36) 


-20 


61 


;o 


29) 


-21 


18 





25) 


-21 


23 





14) 


-21 


37 





11) 





64(0 


05) 





59 


4.80(0.20) 





13(0 


01) 



Note. The Local Group barycenter distances are in Mpc, for Hq = 70 km s~^ Mpc~^ . Absolute magnitudes and g — r colors are corrected for 
internal extinction, and K-corrected including the luminosity distance. The absolute magnitudes are calculated from the SDSS Petrosian magnitudes, 
colors are calculated from the SDSS model colors. Bulge-to-disk ratios, half-light radii, and and disk axis ratios are computed from the i-band GALFIT 
models. Table 2 is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its 
form and content. 

''GALFIT magnitude used instead of SDSS Petrosian, because of poor pipeline estimation of the total observed magnitude. 
''flag-4 galaxies with absolute magnitudes and model colors that do not contain internal extinction corrections and K-corrections. 
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Table 3. Arc-tangent Functional Fit and Rotation Speed 



oUoo name 


I'circ / Km S 


rt 1 arc-sec 


K80 / ^m s 


''end/ Kms 


V 2.2/ Km S 


nag 


J001006. 62-002609. 6 


160.82 


2.79 


129. 27( 2.75) 


138. 57( 2.40) 


125. 90( 2.59) 


3 


J001708.77-005728.9 


126.94 


3.45 


135.45( 6.08) 


125.61( 3.33) 


130. 68( 5.41) 


3 


J001947.33+003526.7 


71.67 


0.84 


85.51(10.59) 


80.07( 5.04) 


85.00(10.42) 


2 


J002025. 78+004934. 9 


123.30 


2.16 


126. 44( 7.67) 


121. 88( 6.43) 


125. 30( 7.51) 


3 


J003112.09-002426.4 


150.47 


1.13 


161.36( 5.46) 


154.90( 3.56) 


159.90( 5.30) 


3 


J003629.24^004342.8 












4 


J003856. 59+001757.3 












4 


J004239.34+001638.7 


224.16 


1.51 


239.59( 3.83) 


240.77( 3.15) 


236. 16( 3.61) 


1 


J004759.81+011209.2 












4 


.1004935.68+010655. 5 


261.47 


15.59 


211.81( 2.24) 


207.32( 1.46) 


200.69( 1.88) 


3 


,1005513.63-010433.8 


217.68 


0.83 


207.71( 5.92) 


205.84( 3.79) 


206.28( 5.44) 


2 


J005516.88-005504.5 












4 


J005603.35-004941.0 












4 


,1005605.99+003742.0 












4 


,)005619.60-000732.7 












4 


J005620. 00-011035. 9 












4 


J005624.61-010834.9 












4 


J005650.61+002047.1 


341.94 


2.66 


310.19( 7.57) 


289.39( 4.66) 


288.43( 5.60) 


3 


J010338.82-002927.9 












4 


J010823.24+000523.1 












4 


J011341. 82-000609.9 


191.40 


1.37 


199.66( 5.84) 


198.76( 3.81) 


197. 77( 5.30) 


1 


J011446.94+003128.9 












4 


J011508.73+003433.6 












4 


J011520. 86+001533. 5 


220.58 


0.33 


212.80(18.88) 


205. 07( 6.21) 


212.83(18.92) 


3 


J011605.60-000053.6 


266.14 


0.81 


255.17(16.61) 


244.61( 7.95) 


252.96(15.14) 


3 


J012000.48-002826.9 


214.53 


0.81 


230.04( 9.13) 


224.99( 5.60) 


228.98( 9.05) 


1 


J012223.78-005230.7 


274.05 


0.85 


290.35( 3.42) 


289.55( 2.51) 


289. 67( 3.09) 


1 


,1012312.37-003828.2 












4 


J012317.00-005421.7 


193.04 


1.93 


176.16( 3.15) 


178.57( 2.37) 


173. 43( 3.01) 


2 


J012340.13+004056.5 


205.76 


1.65 


184.92( 2.37) 


186.13( 1.80) 


180. 95( 2.22) 


1 


J012438.08-000346.4 


196.36 


3.69 


169.58( 3.17) 


170.35( 2.36) 


161. 71( 2.73) 


2 


J012859.57-003342.9 












4 


J013142.14-005559.9 


252.98 


7.29 


227.62( 4.09) 


216.21( 3.29) 


225. 19( 3.92) 


3 


J013600. 15+003948.6 


155.18 


0.89 


180.47( 4.78) 


180.92( 3.20) 


179. 73( 4.74) 


1 


J013752.69+010234.8 


332.39 


2.81 


285.60(39.64) 


242.43( 5.71) 


277.10(34.85) 


3 


J014121. 94+002215. 7 


215.40 


1.64 


202.99( 6.86) 


202. 53( 4.45) 


195.82( 5.73) 


2 


J015550.32+004319.2 












4 


J015746.24^011229.9 


334.50 


0.62 


313.50(13.23) 


310.65( 4.77) 


310.80(11.27) 


2 


J015840.93+003145.2 


178.71 


0.76 


190.29( 2.67) 


189.10( 1.92) 


189.52( 2.63) 


2 


J015852.30+003405.9 












4 


J015946.76-001657.7 


232.66 


1.45 


235.59( 3.86) 


236.24( 3.04) 


230. 93( 3.52) 


1 


J020004.87-001644.8 


186.05 


2.20 


170.36( 2.25) 


175.02( 1.71) 


166. 50( 2.12) 


1 


J020045.13-101451.3 


97.43 


7.35 


83.74(11.32) 


77.50( 5.51) 


75.81( 8.53) 


3 


J020158.20-005951.9 












4 


J020447. 19+005006.3 


108.18 


5.29 


83.77( 4.86) 


78.13( 2.64) 


76.63( 3.64) 


3 


J020540.31-004141.4 


216.65 


0.87 


202.84( 1.80) 


202.81( 1.22) 


201.85( 1.77) 


1 
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Table 3 — Continued 



SDSS name 


V„- „/kms~l 

'Circ/ ^1*^ o 




X^n / IfTM ^ ^ 


' end / ^iii o 






J020621. 28-085218.0 












4 


J020853.01+004712.6 


232.35 


1.74 


239.43( 3.67) 


235. 10( 3.02) 


235. 12( 3.47) 


2 


J021112.60-004933.1 


209.23 


0.78 


223.19( 3.67) 


222.43( 2.65) 


215.97( 3.35) 


1 


J021159.96+011302.8 












4 


J021219.68-004841.4=- 












3 


J021859.65-I-001948.0 


213.66 


2.16 


197.28( 3.75) 


195.41( 2.39) 


191.22( 3.43) 


2 


J021941. 13-001520.4 


94.28 


0.48 


103.09( 4.10) 


103.17( 2.48) 


102.88( 4.08) 


1 


J022026.54-001846.3 












4 


J022606.71-001954.9 


233.91 


2.14 


248.37( 3.89) 


247.41( 2.73) 


245. 95( 3.70) 


2 


J022751.44+003005.5 


252.54 


0.39 


241.67( 2.49) 


242.04( 1.75) 


241. 10( 2.48) 


1 


J022820.86+004114.0 


207.09 


2.51 


219.09( 3.85) 


228.74( 4.54) 


211.34( 3.52) 


1 


J023240.17-I-001536.4 












4 


J023301. 24-1-002515.0 


161.08 


0.69 


181.36( 6.35) 


170.10( 4.12) 


181.65( 6.39) 


1 


J023437.74^005847.3 


269.24 


0.42 


253.44( 3.52) 


250.05( 2.55) 


253.21( 3.51) 


1 


J023522.55-000700.4 












4 


J023610.91-005833.8 


304.27 


3.94 


254.21( 8.57) 


239.53( 5.97) 


240. 75( 6.02) 


3 


J023915.02+002842.5 












4 


J024331. 30+001824.5 


253.07 


1.93 


248.36( 2.48) 


248.93( 2.00) 


246.60( 2.41) 


1 


J024436.70+000831.5 












4 


J024459.89-I-010318.5 


216.47 


3.46 


188.52( 3.48) 


184.25( 2.53) 


180.41( 3.01) 


2 


J024527.56+005451.6 


338.39 


3.84 


313.08( 5.18) 


291. 14( 2.89) 


184.09( 3.19) 


1 


J024558.37-002657.1 












4 


J024727.83-001706.4 


171.50 


1.02 


217.65( 3.75) 


220.52( 3.40) 


215. 77( 3.67) 


1 


J024728.01+003906.9 


216.42 


1.30 


201. 50( 3.95) 


199.21( 2.80) 


196. 17( 3.64) 


2 


J024744.89+002451.0 


185.33 


1.85 


192.08( 3.47) 


190.68( 4.63) 


187.28( 3.27) 


2 


J024850.79-004602.6 


191.65 


0.32 


196.60( 2.86) 


196.75( 1.97) 


196.15( 2.85) 


1 


J024935.34-004325.5 


126.92 


3.68 


103.31( 1.77) 


107.44( 1.77) 


99.89( 1.67) 


2 


J024949.33-004736.6 


230.01 


0.40 


226.87( 2.59) 


227.60( 3.25) 


225. 80( 2.49) 


2 


J025154.59+003953.3 


213.97 


1.93 


206.95( 2.53) 


206.21( 1.84) 


196. 93( 2.22) 


1 


J025338.43+011212.2 












4 


J025627.12+005232.6 


136.03 


3.14 


160.01( 3.89) 


157.02( 2.87) 


157.17( 3.70) 


1 


J025640.44^005614.2 












4 


J030540.35-002406.9 












4 


J032019.21+003005.4 


105.54 


3.89 


91.45(13.88) 


76.41( 4.54) 


83.27(10.19) 


3 


J032251. 25+000822. 7 












4 


J032257.99+000314.9 


227.66 


1.81 


203.71(22.11) 


142.94( 5.37) 


198.88(21.10) 


3 


J033329.46-073308.4 


187.64 


13.25 


113.29( 6.61) 


110.41( 4.97) 


95.73( 3.73) 


3 


J033331. 54+004157.4 


310.93 


4.61 


258.29( 4.96) 


239.79( 3.15) 


252.97( 4.61) 


3 


J033607. 78-003547. 2 


197.80 


2.01 


174.41( 4.08) 


183.96( 3.43) 


174.24( 4.07) 


1 


J035356.55-070422.6 












4 


J080046.85+353146.0 


167.60 


10.17 


106.21(21.90) 


88.71( 5.48) 


89.33(13.12) 


3 


J080126.42+441137.0 


267.69 


2.18 


250.06( 4.74) 


246.41( 2.95) 


246.00( 4.41) 


1 


J080658.75+463346.8 


251.37 


10.54 


179.89( 4.88) 


167.44( 3.47) 


162.05( 3.30) 


2 


J080927.22+435605.6 


202.24 


3.71 


182.20( 3.99) 


183.62( 3.01) 


175.62( 3.67) 


2 


J081313.84+402504.6 












4 


J082949.87+484647.9 


172.61 


5.50 


146.59( 2.22) 


157.94( 1.71) 


138.61( 1.99) 


2 
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Table 3 — Continued 



SDSS name 


VC;,.^ / Icm R ^ 
' Circ / ^3 


Vt 1 3XC-S6C 


Van / Wm S ^ 


K .)/kms~i 

' end / 


V2 2 / km s ^ 


flae 


J082956.27+515824.1 


97.79 


5.33 


86.46( 3.84) 


90.47( 3.54) 


82.26( 3.30) 


1 


J084351. 32+515928. 2^ 












2 


J084408.09+504422.9 


121.74 


6.89 


103.74( 6.04) 


109.18( 3.98) 


96.20( 5.17) 


2 


J084523.21+550707.4 


267.05 


1.59 


261.94( 3.83) 


255. 52( 3.20) 


255. 88( 3.60) 


2 


J085214.30+525323.8 


195.09 


1.26 


213.81( 3.33) 


214.48( 2.51) 


211.84( 3.26) 


1 


J085318.82+570412.5 












4 


J085705.72+514850.7 


118.61 


6.40 


109.71( 5.13) 


100.40( 3.30) 


106.00( 4.68) 


2 


J093814.95+040909.8 












4 


J094949.62+010533.2 


166.84 


3.60 


148.04( 3.42) 


146.78( 2.62) 


140.99( 2.94) 


2 


J095014.42+041155.7 












4 


J095057.66+005147.5 












4 


J095555.07-001 125.0 


117.56 


0.18 


116.53( 6.06) 


116.94( 5.87) 


116.49( 5.97) 


3 


J095743.26+004123.6 


221.69 


2.35 


195.69( 2.89) 


199.62( 2.26) 


188.21( 2.59) 


2 


J100230.82+001826.2 


241.02 


1.17 


260.76( 4.67) 


271. 56( 3.85) 


258. 15( 4.55) 


2 


J100524.08+003238.1 


214.97 


1.87 


207.69( 2.49) 


208.22( 2.35) 


204.30( 2.38) 


1 


J103918.25-001344.8 


101.50 


4.09 


83.97(11.37) 


73.71( 6.11) 


79.14( 9.06) 


3 


J105422.24-002153.9 


330.81 


5.80 


270.95(20.83) 


269.45( 8.20) 


257.60(12.38) 


1 


J111559.50-002059.0 


159.36 


1.82 


160.70( 3.10) 


162.35( 3.25) 


158. 53( 2.92) 


1 


J112346.06-010559.4 


228.74 


3.90 


275.93( 4.62) 


276.66( 4.23) 


258.67( 3.92) 


1 


J115731.83-011510.5 


153.86 


4.50 


137.88( 3.01) 


142.80( 2.43) 


126.54( 2.45) 


1 


J120155.64-010409.3 


210.61 


2.64 


192.54( 5.42) 


193.86( 4.95) 


181.34( 3.68) 


1 


J120430.82+022036.1 


91.10 


3.24 


86.47(18.13) 


75.77( 8.04) 


83.87(16.93) 


2 


J123058.64+513636.2 


75.84 


2.12 


72.25( 6.78) 


70.96( 4.05) 


70.51( 6.38) 


1 


J124428.85-002710.5 


233.38 


6.03 


168.72( 3.65) 


200.76( 4.64) 


159.02( 3.22) 


3 


J124545.20+535702.0 


154.62 


9.63 


100.67( 9.62) 


94.45( 4.35) 


87.77( 6.33) 


3 


J124752.98-011109.0 


189.37 


6.88 


156.93( 2.51) 


165.63( 2.17) 


136.84( 1.80) 


1 


J125554.63-005956.7 












4 


J125715.15-003927.5 


170.06 


2.56 


140.82( 3.47) 


145. 89( 3.04) 


135. 28( 3.11) 


2 


J125854.66-005821.9 


190.55 


0.36 


216.42( 7.73) 


214.91( 6.42) 


204.88( 6.99) 


1 


J130004.72+002710.6 


139.12 


2.14 


114.96( 3.50) 


123.81( 2.84) 


110.06( 3.17) 


3 


J131023.77-010037.9 












4 


J132123.08+002032.7 


229.87 


0.18 


232.18( 2.27) 


232.31( 1.73) 


232.03( 2.26) 


1 


J132903.23-000236.9 


121.19 


1.59 


122.14( 4.36) 


119.45( 3.11) 


119.97( 4.14) 


2 


J133200.24-024255.6 












4 


J133353.35+031735.3 












4 


J133706.48+002204.6 


190.85 


0.20 


225. 78( 4.76) 


226.34( 3.97) 


225. 45( 4.46) 


1 


J133832.17+003225.0 


209.95 


2.17 


201.63( 2.55) 


205.06( 1.91) 


155. 88( 2.40) 


1 


J133839.73+003245.0 


107.80 


3.91 


107.43( 3.99) 


108.12( 2.90) 


104.53( 3.70) 


2 


J134158.85+021844.5 


73.18 


1.95 


72.80(15.83) 


55.70( 5.53) 


70.95(14.57) 


3 


J135102.22-000915.1 


180.52 


0.80 


191.61( 4.95) 


190.83( 2.87) 


190.88( 4.70) 


1 


J135433.67-004635.0 


203.29 


4.12 


162.71( 4.75) 


162. 21( 3.22) 


149.24( 3.67) 


3 


J135946.01+010432.8 


159.76 


0.85 


197.28( 3.73) 


212.71( 3.36) 


195.82( 3.67) 


1 


J140146.46+010924.6 












4 


J140311. 79-005829.4 


175.12 


0.74 


172.25( 6.51) 


170.17( 3.59) 


172.33( 6.54) 


3 


J140342.49-010814.4 












4 


J140407.40+000838.4 












4 
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Table 3 — Continued 



SDSS name 


VC;,.^ / Icm R ^ 
' Circ / ^3 


Vt 1 3XC-S6C 


Van / Wm S ^ 


K .)/kms~l 

' end / ^3 


V2 2 / km s ^ 


flae 


J140452.62-003640.5 


146.75 


2.61 


183.94( 3.83) 


182.62( 3.14) 


180.93( 3.69) 


1 


.1140525.26-004840.9 


244.09 


1.17 


255.57( 4.56) 


259.60( 3.29) 


253. 14( 4.4,3) 


3 


J141026.83-004956.6 


152.92 


0.16 


164.69( 4.57) 


163.58( 3.50) 


164.68( 4.57) 


3 


J141057.69+010208.6 


163.89 


1.01 


154.64( 2. ,36) 


157.13( 1.82) 


151. 72( 2.24) 


2 


J141413.17-005339.8 


262.60 


1.99 


265.06( 2.93) 


263.85( 2.82) 


260.99( 2.75) 


1 


J141628.95-004438.0 


189.62 


1.52 


220.75( 4.04) 


225.96( 3.69) 


217.51( 3.94) 


1 


J142703.51+651154.1 












4 


J142729.65-I-010321.0 


113.14 


2.22 


101.53( 4.11) 


101.62( 2.71) 


98.26( 3.74) 


1 


J143842.98-000027.9 


241.74 


2.85 


189.47( 4.61) 


192.17( 3.03) 


194.20( 4.98) 


3 


J144234.04+005843.9 












4 


J144307.79+010600.0 


192.04 


2.60 


166.32( 3.59) 


165.49( 2.48) 


159.02( 3.13) 


2 


J144418.37+000238.5* 












1 


J144503.29+003137.1 


168.56 


3.42 


158.40( 3.16) 


161.38( 2.19) 


153.51( 2.90) 


2 


J145025.02-011026.5 


271.71 


2.72 


242.75( 3.20) 


241. 79( 2.13) 


234.98( 2.92) 


3 


J150546.86-004253.6 


182.51 


2.44 


173.71( 2.42) 


179.41( 2.21) 


168. 71( 2.24) 


1 


J150645.61+010009.0 












4 


J150739.32-0006.34.0 












4 


J151728.20+010119.4 












4 


J151754.02-I-004943.1 












4 


J152137.86+002757.6 


143.84 


2.00 


134.98( 3.27) 


132.50( 1.87) 


132.54( 3.10) 


1 


J152327.26-010956.0 












4 


J152333.99-002345.3 












4 


J153045.16-002211.5 


211.89 


1.71 


196.18( 2.44) 


199.79( 2.11) 


191.11( 2.29) 


1 


J153914.09+011141.6 












4 


J155452.40+475250.9 


89.30 


3.17 


85.98(25.46) 


66.90( 8.60) 


81.72(22.30) 


3 


J162300. 40+455207.9 


66.57 


2.68 


65.71( 8.74) 


56.03( 5.14) 


62.84( 7.60) 


3 


J170310.47+653417.6 


114.23 


7.10 


81.67( 6.45) 


66.95( 5.65) 


72.07( 3.76) 


3 


J172538.41+592648.3 


127.57 


8.44 


89.73( 4.57) 


94.46( 2.73) 


80.06( 3.38) 


1 


J172747.29+541225.4 












4 


J203523.80-061437.9 


311.82 


15.89 


258.39( 3.13) 


244.96( 1.78) 


249. 69( 2.78) 


3 


J203755.29-061800.2 


245.29 


3.56 


224.69( 8.38) 


225.70( 4.47) 


221.98( 8.15) 


1 


J203907.00+003316.2 


172.82 


0.39 


161.45( 3.87) 


159.78( 2.75) 


160.31( 3.77) 


3 


J204012.01-055546.4 












4 


J204036.79-044133.8 


82.40 


2.22 


66.84( 6.00) 


63.06( 3.89) 


63.72( 5.09) 


3 


J204256.27-065126.1 


135.43 


2.01 


118.86( 5.66) 


115.50( 3.26) 


114.40( 5.05) 


1 


J204515.04-053721.8'' 












1 


J204805.63+000407.4^ 












2 


J204813.52-054741.7'' 












2 


J204913.40+001931.0 


236.91 


5.04 


179.69( 9.34) 


174.42( 5.69) 


140.32( 4.63) 


2 


J205103.70+000825.5 


179.13 


8.35 


125.08(13.75) 


110.11( 4.42) 


110.24( 9.94) 


3 


J205212.77-054754.0 


153.04 


7.71 


124.05( 2.48) 


125.92( 2.73) 


115. 58( 2.08) 


2 


J205235.43-054240.0 


322.05 


6.59 


297.73( 6.17) 


280.28( 3.82) 


288. 67( 5.,34) 


1 


J205307.50-002407.0 


121.00 


3.60 


116.05( 5.06) 


112. 55( 3.08) 


110.66( 4.38) 


3 


J205404.34+004638.6 


222.30 


2.65 


215.31( 3.57) 


215.61( 2.42) 


210.79( 3.31) 


1 


J205432.43+000231.6 


179.32 


1.18 


220.23( 5.06) 


215.19( 3.91) 


218.18( 4.95) 


2 


J205528.79-071951.3 


80.13 


1.57 


76.04( 8.76) 


70.99( 4.70) 


73.98( 8.06) 


3 
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Table 3 — Continued 



SDSS name 


VC;,.^ / Icm R ^ 
' Circ / ^3 


Vt 1 3XC-S6C 


Van / Wm S ^ 


K j/kms~l 

' end / ^3 


V2 2 / km s ^ 


flae 


J205532.62+000635.6 


142.62 


5.38 


105.83( 3.71) 


107.33( 3.23) 


95.67( 2.79) 


2 


.1205755.40+000711.5 


1,35.56 


1.18 


130.13( 3.50) 


133.80( 2.81) 


128. 41( 3.37) 


2 


.1210039.64-001236.6 


126.82 


4.34 


105.43( 3.49) 


102.06( 2.96) 


99.77( 2.90) 


3 


.1210551.33+091237.4 












4 


J210633.54+104504.1 


110.00 


2.44 


108.67( 4.31) 


105.13( 2.94) 


105.05( 3.86) 


2 


J211222.69-074913.9 


85.96 


3.54 


73.35( 4.32) 


68.39( 3.19) 


72.19( 4.08) 


2 


J211343.93+003428.7 


317.25 


3.35 


279.42( 4.85) 


276.96( 4.81) 


270.41( 3.24) 


2 


J211439.91-075806.9 


119.26 


2.36 


111.01( 4.02) 


108.61( 2.63) 


107.71( 3.66) 


1 


J211450.23-072743.3 


182.75 


1.57 


213.79( 4.77) 


211. 17( 3.78) 


211. 47( 4.67) 


1 


J211522.10-074605.0 


68.94 


1.77 


58.69( 5.97) 


56.82( 4.41) 


54.43( 4.66) 


2 


J211601. 74-074423.1 












4 


J211627.63-004935.3 


285.69 


9.94 


286.33(14.63) 


236.98( 4.63) 


269.32(12.28) 


3 


J211816.06-073507.8 


333.30 


10.46 


241.95( 7.34) 


237.06( 4.57) 


222.01( 5.48) 


3 


J211900.28-005750.0 


144.06 


0.98 


153.46( 3.82) 


158.74( 3.20) 


153.91( 3.85) 


1 


J213014.23-080401.1 


191.32 


1.03 


191. 74( 3.04) 


190.32( 1.94) 


190.64( 2.87) 


1 


J213058.11-070507.9 


159.43 


0.85 


154.61( 6.74) 


145.92( 2.85) 


153. 30( 6.60) 


1 


J213703.87-073518.0 


141.73 


1.68 


136.19( 3.79) 


138.14( 2.73) 


131.44( 3.43) 


1 


J213811. 68+121139.1 


113.68 


2.31 


101. 74( 3.87) 


99.43( 3.00) 


96.98( 3.32) 


2 


J214413.21+114936.4'' 












2 


J215156.74+121411.3 


145.64 


2.90 


122.23( 5.87) 


115.99( 3.28) 


116.42( 5.05) 


3 


J215226.03-081024.9 


203.56 


0.80 


190.58( 4.41) 


184.59( 2.81) 


188.88( 4.28) 


2 


J215247.62+122942.8 


129.53 


2.84 


124.49(23.03) 


91.27( 6.99) 


120.95(20.26) 


3 


J215326.90+002218.0 


131.27 


3.19 


135.90( 3.95) 


135. 17( 2.69) 


132. 13( 3.65) 


1 


J215421. 67-075605.7 


163.89 


2.09 


135.73( 2.69) 


148.60( 2.55) 


126.27( 2.51) 


2 


J215435.36+125140.0* 












3 


J215652.70+121857.5 


161.72 


1.19 


180.06( 3.55) 


179.94( 2.88) 


178.49( 3.46) 


1 


J215837.69+114913.2 


188.98 


2.40 


206.42(10.87) 


195.51( 5.26) 


200.25(10.09) 


2 


J230857.95-002302.9 












4 


J232021. 17-001819.2 


148.71 


6.46 


128.90(81.92) 


87.12( 6.93) 


117.27(71.37) 


3 


J232238.68-005903.7 


161.82 


3.38 


168.56( 3.48) 


164.61( 2.94) 


162.72( 3.15) 


2 


J232338.75+011334.8 












4 


J232613.88+010828.2 


183.74 


3.08 


166.23( 3.49) 


165.85( 2.72) 


158.25( 2.99) 


2 


J232631. 10+005013.5 


182.75 


2.55 


165.44( 2.89) 


175.59( 2.68) 


157.95( 2.59) 


2 


J233032.04+000451.7 












4 


J233152.99-004934.4 


121.83 


9.22 


88.72( 4.27) 


88.10( 2.55) 


80.27( 3.21) 


2 


J233259.3,3+004318.8 


113.23 


6.92 


86.72(11.26) 


87.53( 6.96) 


78.70( 8.09) 


3 


J233908.64+002748.3 












4 


J234107.77-010215.2 


222.56 


0.94 


217.13( 2.32) 


226.81( 1.94) 


211.56( 2.25) 


1 


J234125.92-003943.9 












4 


J234215.64^003745.9 


180.17 


2.25 


157.19( 4.01) 


158.12( 3.21) 


151. 32( 3.58) 


3 


J234328.26-000148.6 


226.11 


2.51 


194.95( 3.73) 


196.62( 2.48) 


187.09( 3.29) 


2 


J234504.86-001615.1 


131.65 


3.84 


117.40( 2.72) 


115.98( 1.90) 


113.63( 2.45) 


2 


J235106.25+010324.0 


253.79 


10.21 


223.41( 6.29) 


214.22( 3.70) 


212.25( 5.02) 


1 


J235410.09+002258.3 












4 


J235603.89-000958.6 


114.18 


3.26 


102.95( 2.83) 


104.38( 2.21) 


99.18( 2.55) 


2 


J235607.82+003258.1 


179.72 


0.90 


233.57( 6.84) 


235.00( 5.38) 


232.03( 6.74) 


2 
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Table 3 — Continued 



SDSS name 


Vcirc/kms ^ 


rt 1 arc-sec Vso / km s ^ 


Vend/ km s 1 


V2.2/kms ^ 


flag 


J235624.68-001739.6 


112.94 


2.94 96.80( 9.61) 


83.80( 4.17) 


90.78( 7.80) 


3 


,1235631.52+003300.1 










4 


J235656.66-005912.3 










4 


J235743.74+003918.6 


192.70 


2.92 168.35( 2.27) 


173.91( 1.71) 


161.25( 2.03) 


2 



Note. — Parameters of arc-tangent functional fits after 10 kms~^ has been added to the velocity measure- 
ment errors on each data point. The velocities have been corrected for inclination and cosmological broadening. 
Uncertainties are listed in parenthesis, flag-4 galaxies were observed but did not have enough Ha: to determine a 
rotation curve. Tabic 3 is published in its entirety in the electronic edition of the Astronomical Journal. A portion 
is shown here for guidance regarding its form and content. 

''Rotation curve data not used because of a poorly determined SDSS pipeline isophotal position angle. 

''Rotation curve data not used because of a poorly determined SDSS pipeline isophotal axis ratio. 
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Table 4. TF Parameters 



y 


X 


xo 


a(±) 


a'(±) 


6(±) 


6'(±) 


a{±) 


a'{±) 












Vso flag 


:-l,2,3 N =162 










logVso 


2.220 


-5.476(0.226) 




-20.686(0.041) 




0.457(0.034) 




logVso 




-20.607 


-0.136(0.006) 


-7.35(0.32) 


2.209(0.006) 


-20.687(0.049) 


0.073(0.006) 


0.536(0.044) 


Mr 


logVgo 


2.220 


-5.962(0.204) 




-21.182(0.035) 




0.418(0.036) 




logVso 


Mr 


-21.107 


-0.135(0.006) 


-7.41(0.32) 


2.210(0.006) 


-21.181(0.044) 


0.063(0.005) 


0.466(0.037) 


Mi 


logVso 


2.220 


-6.321(0.220) 




-21.390(0.035) 




0.423(0.035) 




logVso 


Mi 


-21.327 


-0.130(0.005) 


-7.69(0.29) 


2.212(0.005) 


-21.388(0.038) 


0.061(0.005) 


0.469(0.038) 


M^ 


logVso 


2.232 


-6.590(0.238) 




-21.518(0.034) 




0.415(0.038) 




logVso 


M^ 


-21.400 


-0.127(0.005) 


-7.87(0.31) 


2.216(0.005) 


-21.525(0.39) 


0.058(0.006) 


0.457(0.047) 



Note. — Bivciriate relations are fit with the model y=a{x-xo)+b with a Gaussian intrinsic scatter of y at fixed x, with dispersion 
(7. Errors listed for a, b, and a are computed from 100 bootstrap trials, and the value of ceo is chosen so that errors in a and b 
are uncorrelated. a' , b' , and a' are the forward equivalent of the inverse relations. 
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Table 5. Bivariate Relations and Residual Correlations 



y 


X 


xo 


a(±) 


b(±) 


a(±) 


g-r 


logVgo 


2.220 


0.576(0.055) 


0.552(0.0079) 


0.108(0.006) 


r (V-80) 


AMg (Vso) 


0.00 


0.067(0.018) 


0.000(0.0075) 


0.102(0.0067) 


g-r 


Mg 


-20.607 


-0.063(0.018) 


0.544(0.0069) 


0.128(0.0062) 


Ag-r {Mg) 


AlogVgo [Mg) 


0.00 


0.806(0.122) 


0.007(0.007) 


0.091(0.0069) 


Ag-r (Vso) 


AM, (Vso) 


0.00 


-0.003(0.021) 


0.000(0.008) 


0.108(0.006) 


g — r 


M, 


-21.327 


-0.075(0.008) 


0.547(0.009) 


0.115(0.008) 


Ag-r(M,) 


AlogVso (M.) 


0.00 


0.534(0.151) 


0.000(0.008) 


0.108(0.007) 




logVgo 


2.220 


0.758(0.074) 


0.755(0.0115) 


0.169(0.010) 


AlogRi 


AM, 


0.00 


-0.130(0.036) 


0.000(0.011) 


0.155(0.0097) 


log/?. 


M, 


-21.327 


-0.121(0.011) 


0.749(0.011) 


0.157(0.009) 


Alogi?, {AU) 


AlogVso (M,) 


0.00 


-0.046(0.195) 


0.00(0.011) 


0.157(0.009) 



Note. — 



For modeling procedure, see note to Table |4] 
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Table 6. Alternative TF Fits 



y 


X xo o-i^) 


6(±) 


a{±) 






Vso N = 162 Weighted by fia 


ig-1,2,3 Mr distribution 







logVgo 


2.22 


-5.492(0.272) 


-20.671(0.042) 


0.457(0.041) 


lug VSO 


M 

IVlg 


-20.607 


-0.136(0.007) 


2.204(0.007) 


0.073(0.007) 


Mr 


log Vso 


2.22 


-5.977(0.268) 


-21.166(0.0392) 


0.418(0.037) 


1 C\rT\/ci r\ 


Mr 


-21.107 


-0.135(0.006) 


2.208(0.00610) 


0.064(0.007) 




logV^o 


2.22 


-6.334(0.286) 


-21.347(0.039) 


0.423(0.038) 


logV^o 


Mi 


-21.327 


-0.130(0.005) 


2.211(0.006) 


0.061(0.006) 


M 


log vso 


2.22 


-6.607(0.260) 


-21.496(0.044) 


0.415(0.041) 


logV^o 


IVlz 


-21.400 


-0.127(0.005) 


2.216(0.007) 


0.058(0.007) 




Vso flaj 


;-l,2,3 N = 162; No Internal Extinction Correction 




log Vso 


2.220 


-5.190(0.200) 


-20.174(0.039) 


0.454(0.029) 


''80 


M 


-20.607 


-0.145(0.006) 


2.284(0.007) 


0.075(0.005) 


Mr 




2.220 


-5.719(0.189) 


-20.775(0.030) 


0.418(0.033) 


log VQO 


Mr 


-21.107 


-0.141(0.006) 


2.268(0.006) 


0.065(0.005) 


M' 


log K80 


2.220 


-6.166(0.217) 


-21.092(0.033) 


0.413(0.033) 






-21.237 


-0.134(0.005) 


2.253(0.005) 


0.060(0.006) 


Mz 


log V go 


2.232 


-6.438(0.231) 


-21.231(0.029) 


0.408(0.039) 


logVso 


Mz 


-21.400 


-0.130(0.005) 


2.243(0.005) 


0.058(0.005) 






V2. 


.2, fiag-1,2,3 N = 


162 




Mg 


logV2.2 


2.220 


-5.420(0.213) 


-20.777(0.038) 


0.438(0.030) 


logV2.2 


Mg 


-20.607 


-0.142(0.006) 


2.191(0.006) 


0.072(0.006) 


Mr 


logV2.2 


2.220 


-5.907(0.196) 


-21.282(0.032) 


0.392(0.032) 


logV2.2 


Mr 


-21.107 


-0.140(0.005) 


2.192(0.006) 


0.061(0.005) 


Mi 


logV2.2 


2.220 


-6.263(0.198) 


-21.498(0.033) 


0.391(0.032) 


logV2.2 


M, 


-21.327 


-0.135(0.005) 


2.195(0.005) 


0.058(0.005) 




logV2.2 


2.220 


-6.507(0.219) 


-21.548(0.032) 


0.391(0.033) 


logV2.2 


Mz 


-21.400 


-0.131(0.004) 


2.199(0.005) 


0.056(0.005) 






Vend flag-1,2,3 N = 


162 




Mg 


logVcnd 


2.233 


-5.309(0.207) 


-20.789(0.039) 


0.461(0.035) 


logKnd 


Mg 


-20.363 


-0.144(0.006) 


2.164(0.007) 


0.077(0.006) 


Mr 


logVend 


2.207 


-5.773(0.189) 


-21.142(0.036) 


0.432(0.036) 


logVend 


Mr 


-20.873 


-0.141(0.006) 


2.169(0.006) 


0.068(0.005) 


Mi 


logVend 


2.214 


-6.109(0.219) 


-21.392(0.035) 


0.442(0.035) 


logVend 


Mi 


-21.122 


-0.136(0.006) 


2.176(0.006) 


0.066(0.005) 
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Table 6 — Continued 



y 


X 


XQ 


a(±) 


b{±) 


<7(±) 




logVcnd 


2.212 


-6.362(0.224) 


-21.428(0.037) 


0.439(0.034) 






-21.210 


-0.133(0.005) 


2.182(0.005) 


0.063(0.006) 








Vso, flag-l,2N = 


114 




Ms 


logVso 


2.220 


-5.157(0.326) 


-20.776(0.050) 


0.431(0.039) 


logVso 




-20.607 


-0.137(0.008) 


2.202(0.008) 


0.071(0.007) 


Mr 


logVso 


2.220 


-5.599(0.299) 


-21.259(0.044) 


0.396(0.038) 


logVso 


Mr 


-21.107 


-0.137(0.008) 


2.203(0.008) 


0.063(0.007) 




logVso 


2.220 


-6.024(0.323) 


-21.455(0.043) 


0.401(0.039) 


logVso 


M, 


-21.327 


-0.132(0.007) 


2.208(0.006) 


0.059(0.006) 




logVso 


2.232 


-6.266(0.382) 


-21.566(0.048) 


0.391(0.042) 


logVso 


M^ 


-21.400 


-0.130(0.008) 


2.213(0.006) 


0.056(0.007) 



Note. — Same as Table U but with different weigliting of data points, internal 
extinction correction, velocity definitions, or sample selection, as indicated. 
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Table 7. Tully-Fisher Relations in the Literature 



study 


data 


bias corrected 


slope 


scatter 


PTSS^' 


i?-HI 


no 


-7.64(0.29) 


< 0.25 mag 


PT92'' 


R-HI 


no 


-8.23 


small 




r-Ho 


yes 


-8.195(0.263) 


0.3 


V01<^ 


R-HI 


no 


-9.6(0.4) 


0.15 


K02° 


i?-band-Ha 


no 


-9.70(0.27) 


0.42 


PT88'' 


7-HI 


no 


-7.85(0.29) 


< 0.25 mag 




7-HI 


no 


-8.72 


small 


W95f 


7-HI 


no 


-7.865(0.162) 


0.392 


WQSf 


7-HI 


yes 


-6.925(0.158) 


0.360 




7-Ha/HI 


yes 


-6.80(0.08) 


0.3 


G97S 


7-HI 


no 


-7.16 


0.2-0.3 mag 


G97g 


7-HI 


yes 


-7.52 


0.2-0.3 mag 


C97h 


r-H« 


no 


-6.36(0.22) 


0.46 mag 


C97' 


7-Ha/HI 


no 


-6.64(0.13) 


0.56 mag 



Note. — Tully-Fisher relations from different samples. The intercepts 
are not listed due to different assumptions for 77o, and different values for 
the assumed internal extinction corrections. The scatter is the intrinsic 
scatter. 

^Fierce fc TuUvl l IlQSSi) study of spirals in Virgo and Ursa Major, having 
no morphological peculiarities, depth of the cluster included in analysis. 

^Pierce &: TuUvl l ll992h study of galaxies with Cephied/RR-Lyrae, and 
PNe distances. Sample contains 15 galaxies. 

' IWillick et"al] l ll996l) analysis of the CF sample, with cz > 3000 km s"! 
in the Local Group frame, with the total scatter. 

■^VOl forward TF relation for the RC/FD Sample (excluding NGC 
3992), using Vfiat, and the most likely intrinsic scatter (depth of clus- 
ter assumed to be 0.17 mag). 

'^K02 inverse TF relation for galaxies with A/jj < —18, weighted inverse 
slope. 

IWillick et al] lll995l ) TF relations for spiral galaxies in clusters. 
e lGiovanelli et al.l lll997l) cluster TF relation for 782 spirals. 
''C97 analysis of CF data with optical rotation curves. 
'C97 analysis of MAT sample having optical and HI data. 



